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CHAPTER  I 


INTRODUCTION 

All  communication  channels  suffer  degradation  from  nonlinearities. 
Analysis  and  correction  of  the  nonlinear  distortion  thus  produced  is 
seldom  feasible  with  just  linear  methods,  since  they  fail  to  capture  all 
of  the  phenomena  involved.  Because  of  this,  the  Nonlinear  Transfer 
Function  (NLTF)  approach,  based  on  the  Volterra  theory  of  nonlinear 
systems  [1],  has  recently  been  applied  for  weakly  nonlinear  systems  [2]. 
A  weakly  nonlinear  system  is  one  for  which  the  response  is  dominated  by 
the  linear  contributions  for  (amplitude)  ranges  of  the  input  which  are 
encountered  in  practice.  Such  systems  are  also  sometimes  referred  to  as 
quasi-linear  systems. 

In  communication  systems  many  sources  for  nonlinear  distortion 
exist  from  nonlinear  discrete  devices  in  amplifiers,  metal-to-metal 
oxide  junctions  in  aircraft  shells,  antenna  structures,  etc.  The 
susceptibility  to  Interference  can  be  especially  pronounced  when  the 
transmitter-receiver  pair  is  located  in  an  electronically  dense  command 
platform  having  a  profusion  of  RF  emitters  and  receptors.  For  example, 
when  multiple  carriers  are  amplified  simultaneously  by  one  transmitter, 
Intermodulatlon  (IM)  products  are  generated  due  to  the  nonlinearities  in 
the  power  amplifier  (TWT  or  Klystron).  Similarly,  a  strong  (locally) 
transmitted  signal  leaking  into  a  receptor  can,  when  processed 
simultaneously  with  a  weak  but  desired  received  carrier  by  a  nonlinear 


V  ■  V  "L^  ''VT  H,"  •  •'•'I*’*  ■’1*  I!  *1.'  W  ^  *  'J  •fl  ■'■. 


element,  produce  intermodulation  effects. 

This  report  introduces  a  practical  method  for  compensating  the  un¬ 
desirable  nonlinear  effects  produced  within  such  a  communication  link. 
Using  the  known  characteristics  of  the  channel  and  an  appropriately 
selected  post-compensator  structure,  optimum  parameter  values  are  found 
for  the  compensator.  The  resulting  design  achieves  significant  reduction 
of  nonlinear  effects. 

As  stated  above,  we  cast  the  compensation  problem  in  terms  of  a 
Mean-Square  (MS)  intermodulation  criterion.  This  MS-IM  criterion, 
formulated  for  a  frequency  band  of  interest,  is  minimized  by  the  design 
program  VC0MP3.  This  program  yields  the  optimum  compensator  parameters. 
High  reliability  in  the  minimization  process  is  achieved  by  use  of  a 
powerful  optimization  package  NL2SN0  [M].  Clever  parametrization  of  the 
block  Transfer  Functions  (TF)  guarantees  the  stability  of  the  complete 
compensator  over  the  entire  parameter  space.  Examples  provided  in  the 
report  demonstrate  both  the  simplicity  and  the  high  degree  of 
effectiveness  achievable  through  this  new  methodology. 
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CHAPTER  II 


REVIEW  OF  VOLTERRA-BASED  SYSTEMS 


In  the  analysis  of  wide  band  amplifiers,  it  is  often  assumed  that 
the  output  depends  only  on  the  input  signal  applied  at  the  same  instant 
of  time.  The  input/output  relation  can  thus  be  expressed  with  a  power 
series  expansion  as  follows: 

y(t)  -  a^ x(t)+a2x^(t)+a^x^(t)+  •••  (1) 

where  x(t)  and  y(t)  denote  the  input  and  output  signals,  respectively, 
and  the  coefficients  a^  are  time- independent  constants.  In  general, 

however,  the  output  y(t)  is  also  dependent  on  the  past  input  signal.  An 
output  expression  that  Includes  the  contribution  of  the  previous  input 
variations  is  the  sum  of  multidimensional  convolution  integrals 


y„(t)  =//.../  h^(x,,T2,--'.T„)x(t-T,)x(t-T2) 

—  00  —00 

...x(t-T^)dT,dT2  •••  dx^  (3) 

and  h^  is  a  real-valued  symmetric  function  of  n  real  variables. 

Expression  (2)  is  usually  referred  to  as  the  Volterra  series. 

This  representation  shows  that  a  nonlinear  system  may  be  regarded  as  the 
combination  of  a  linear  and  a  number  of  higher  order  nonlinear 


subsystems.  Each  of  these  subsystems  is  characterized  by  a  n- 


dimensional  impulse  response  ’  ^2’ *  *  *  *  '  also  called  the  nth  order 

Volterra  kernel  [1].  For  a  physically  realizable  system,  h^  is  causal, 

i.e.,  it  has  the  value  zero  whenever  any  of  its  arguments  is  negative. 
Also,  these  kernels  must  be  absolutely  summable  for  a  stable  system. 

The  nth  order  transfer  function  is  defined  as  the  n-fold  Laplace 
transform  of  h^,  i.e.. 


•'"n^ 


;  / 


/  h^(T,.T2. 


•.Tn)e 


-(SiT^* 


+  3  T  ) 

n  n 


•  dT,dT«»**dT  (4) 
I  2  n 


In  particular,  we  shall  call  H^(3^)  the  linear  transfer  function. 

In  Fig.  1  a  block  diagram  is  shown  representing  the  multi-order 
responses  summed  to  produce  the  final  system  response. 


1  Fig.  1 .  Model  of  a  weakly  nonlinear  channel  using  the 

I 


Volterra  series  representation. 


VOLTERRA  REPRESENTATION  OF  THE  QUADRATIC  SYSTEM 


In  this  section  we  will  discuss  in  detail  the  second-order  (i.e., 

quadratic)  Volterra  subsystems.  Fig.  2  shows  an  example  of  a  quadratic 

subsystem  where  H  (s),  H^(s)  and  H  (s)  are  linear  transfer  functions, 
a  b  c 

This  particular  example  represents  the  most  basic  second-order 
subsystem  and  can  be  shown  [2]  to  be  characterized  by  the  quadratic 
transfer  function 

H-(s, .s_)  =  H^(3, )H.  (s„)H^(s  +s„)  (5) 

212  a1b2c12 

Note  that  in  general 

*^2^  ^1  ’  ^2^  *  ^2’  ^1  ^ 

where  h^  denotes  the  unsymmetric  form  of  the  kernel.  The  unsymmetric 
kernel  can  be  symmetrized  by  defining  a  new  kernel  as 

^2  ^  ^  \  ^  2  ^  ^2  ^  ^1  ^  ^2  ^  *  ^2  ^  ^2  ^  ^1  ^  ^ 

Since  it  is  not  possible  to  measure  h2(T^.T2)  only  the  system  input 

and  output  [1],  it  makes  sense  to  consider  only  symmetric  kernels  in  our 

analysis  without  any  loss  of  generality.  In  particular,  we  shall  assume 

that  H  (s)  »H^(3)  for  the  second-order  subsystem  shown  in  Fig.  2. 
a  b 

Correspondingly  the  time-domain  Volterra  kernel  can  be  shown  to  be  [2] 

00 

h  (a.  ,0^)  ==  /  h  (T)h  (o  “T)h  (o  -T)dT  (6) 

2  12  c  a  1  D  2 

—  00 

Associated  Two-dimensional  Response 

The  expression 

00 

^(2)  ^''l ’''2^  =  /  /  h^it^  ,T2)x(t^-Tpx(t2-T2)clT^dt2  (7) 


is  known  as  the  associated  two-dimensional  response.  Its  significance 


arises  from  the  fact  that  if  we  set  t^»t2=t  then, 

y(2)(t,t)  -  y2(t) 


(8) 


Fig.  2.  Basic  second-order  Volterra  system. 


That  is,  contains  the  response  y^i^)  along  the  ^45°  line  in 

the  t^-t^  plane.  Further,  this  associated  response  can  be  computed 

through  linear  operations  on  the  input  signal.  Indeed,  the  two- 
dimensional  Fourier  transform  of  y(2) *''2^ 

,a)2)X(a)^  )X((i)2)  (9) 

where  X((d)  is  the  Fourier  transform  of  the  input  signal.  Taking  the 
two-dimensional  inverse  Fourier  transform  yields 


(2it) 


(10) 


Hence,  by  setting  t^-t^-t,  we  have 


dip2 


y2(t)  -  I  I  )X(*2) 

(2ir)  “« 


di)»^  (11) 


Further,  by  taking  the  Fourier  transform  of  (11)  we  obtain 


^2^“^  “  ~2Ti  ^  ^  H2(^^  ,\|<2)X(il»^  )X(i|)2)6((»)-i|»^ -tti2)d4'^<li|)2 


In  view  of  (9)  we  can  also  write 


“  ~2^  /  X(2)^'>’l*“"*l)<^'^1  ^’3) 

“■  OB 

In  other  words,  ¥2(0))  is  the  integral  of  ^^2)  ^*•'1  >’*'2^  along  the  line 
+i|)2:  alternatively,  Y^iu)  is  the  profile  of  ^^2)  •’•’2^  along  the 

line  4'^=4'2* 

Let  us  now  consider  two  simple  examples  in  order  to  better 
appreciate  the  intermediate  use  of  the  associated  response  for  finding 


vy..;.v; 


Cv.'v'.v-.-. 

'.'-V-V-'.- 


*•  •  •  V  • 

•  »  '  •  • 


the  nonlinear  system  response. 


Example  1 


Let  us  compute  the  quadratic  system  response  to  a  single  complex 
sinusoid. 


x(t)  ■»  Ae 

Translated  into  the  frequency  domain  the  input  is  written  as 


X( u)  »  2irA6( u)~n^  ) 


k  k  »  «  •  »  « 


Y^2)  (“i  »‘^2^  ”  H2(u)^  ,(D2)'5(u)^~£J^  )6(a)2"np 


•  tr  S  f  H2('|/^  ^'^^“''’•'2”’*’l  ^^’*'2^'*'l 

—  00 

-  7  - 


-.'.A'. 


yy\'< 


A  -*vV-  -’-V*'  .• 

I*  «**  /  \ 


.* 


Then 


—  0» 

-  A^H2(a,.a^)6( (1.-20^) 

j2n  t 

y^ct)  -  ,Q,  )e  (16) 

We  note  that  the  output  frequency  of  the  quadratic  nonlinearity  is  twice 
the  input  sinusoids's  frequency. 

Example  2 

Suppose  now  that  the  input  to  the  quadratic  system  consists  of  the  sum 
of  two  sinusoids,  given  by 

jfl.  t  jfi  t 

x(t)  -  A^e  ^  (17) 

X(a))  =■  2T<[k^6(.ui-a^)->-k^6(.ui-a^)] 

The  corresponding  output  is 


—  09 

•  [a^  -i|)2)di))^  dij)^ 

p  j2n  t  j(n  +n2)t 

y2(t)  =  A^H^Cn, ,n^)e  ♦2A^ A^H^C , 02)6 

P  J2a  t 

+A2H2(n2.«2^® 

Thus  an  input  consisting  of  the  sum  of  two  sinusoids  results  in  an 
output  containing  three  frequencies. 

The  results  of  the  above  examples  can  be  generalized  to  the  case  of 
an  input  consisting  of  the  sum  of  M  sinusoids,  expressed  by 


M  jn,t 
x(t)  -  I  A  e 


vv** 


j-  > 


(19) 


The  quadratic  response  to  this  input  can  be  shown  to  be 


MM  2!  A, A, 

y2(t)  -  I  I  - - 

1  S  k  S  i  tn^  •  •  •  n^j! 


H2(Q^,fl^)e 


j(Q^+Oj^)t 


where 


”l  '^”2'^  *  ’  *  ^ 

and  m  is  the  number  of  occurrences  of  the  £th  frequency  (as  arguments  of 


SECOND-ORDER  VOLTERRA  SYSTEM  RESPONSE  TO  THE  SUM  OF  M  REAL 

SINUSOIDS 


Multi-channel  communication  systems  must  necessarily  employ 
multiple  carrier  waveforms.  In  a  linear  environment  no  new  frequencies 
are  produced  which  were  not  present  in  the  input.  Nonlinearities,  even 
of  a  mild  nature,  will  generate  new  terms.  We  derive  a  closed  form 
expression  for  the  response  of  a  second-order  system  to  an  input 
consisting  of  the  sum  of  M  sinusoids.  Given  the  second-order  symmetric 
nonlinear  system  of  Fig.  3,  let  the  input  be  the  sum  of  M  cosine 


waveforms 


x(t)  -  A^  cos(Q^  t)+A2C0s(n2t)+  •••  +Aj^cos(Q|^t) , 


The  output  at  v^  is  denoted  by 


V  (t)  ■  V  +v  +  . • .  +v 

a  a^ 


AjHa(n^)  |cos(flit  ♦Sjj  )+A2|Ha(n2)  jcosCn^t  )*  ••• 

^  •  ^22) 

M 


where  e„  is  defined  as  Phase{H  (£2,)} 

a  1 
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To  proceed  further  with  the  analysis,  let  us  compute  the  response 


at  w, 


W(t)  “  V  ‘V 

a  a 


M  M 


■I  I  V  V 
a.  a, 

i»1  k-1  1  k 


M  M 

I  1 


■■  I  I  -  V  V 

1  <  i  S  k  •••  n^j!  ^k 


Fig,  3.  Basic  second-order  nonlinear  system. 


where 


.  2!A,A, 

A  i  k 


"  *"2’  *  *  *  **Vl' 


and  0-  =  Phase{H  (a,)} 
£Jj  a  i 


Equation  (23)  may  be  simplified  by  using  the  trigonometric  identity 


cos(  t)cos(  ^2^)  -  \  [cos(5^+52^*^  +cos(5^-52)^] 
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•  '  *•  , 
C.'-V*/* 

.V ‘.Vr 


KvI-Iv/J 


•  ,  -  V  a  * 


*  ** ***  ****  **■*  ^'*i> 


-Ik 

fv 

!i 


Thus,  we  have 


M  M 


““’  '  ^  l^s  k  S  i‘’lKl"a‘“l>WH°"'‘W’- 


'“=‘<“l-“k>t  *®0  -»0 
1  k 


The  output  is  then  written  as 


M  M 


'''*  '  \  .1  .1  .  ■'lkl"a<‘‘l>WI 


1  ^  k  S  i 


•lH<,(ai-VI“=«V“k>‘  ‘V-VVa.” 

i  k  k  1 


To  illustrate  the  use  of  this  consider  an  input  consisting  of  the  sum  of 


two  sinusoids, 


K(t)  -  C03(ft^ t)+A2COs(«2t) 


The  output  is  then 


y(t)  -  -  Af|H^(n^)p|H^(2a^)|[co3(2n,t  *620^ )*i  ] 


^A,A2|H^(n^)H^(n2)H^(n,*Q2)l[cos((R,^a2)^ 


^cos((«^-«2)t 


^  2  *2l“a^”2^1  |Hc(2n2)|[co3(2fi2t  +  ]  (26) 


where  41^^  «  Phase{H^(  fl j^) ) 


Since  H  (a  , a  )  «  H  (a, )H  (a  )H  (a  +a„),  the  above  expression  simplifies 

ci  \  i-  ^  I  Si  ei  C  i  <L 


y(t)  -  ^  A^[|H2(a,,a,)lco3(2a,t  )+Hp(a,  ,-a, )  ] 


I  v»  '  /  '  1*2  '  t 


.  JK. 


mm 


*•  • 
.  •  .  •  -  •  %  • 


**• 


rZ-’S:*’.’ 


.  - . .  .-.'.v 


+A^ |coa((ft^  *“2^ 

♦|H2(n^ ,-ft2)|cos((a^-ft2)t  +a^)] 

*  ^  A2[|H2(n2.«2Hcos(2Q2t  )  +  | H^C . -«2 H  ^  (27) 

where 

=  PhaselH^Cf.j  ,B^  ) } 

=  Phase{H2(B^ .B^) 1 
=  Pha3e{H2(B^ ,-fl2) ) 

=  Phase{H^(B2,B2)} 

THIRD-ORDER  NONLINEAR  RESPONSE  TO  A  SUM  OF  MULTIPLE  COSINE  CARRIERS 


The  cubic  nonlinear  response  to  a  sum  of  M  cosine  carrier  waveforms 
is  computed  as  above.  Given  the  nonlinear  system  of  Fig.  4  the  response 
to  an  input  of  M  cosines  is 


y(t)  =  i  ^  I  I  r  |Q  1 
i  i  k  <  I 


M  M  M 


(28) 


ikJ, 


where 


A 


^ik*, 

**  '*•  -* 

.*  vV. 

'^ikt 

'«i 

-  Phase{H  (ft,)} 
a  1 

:  m 

®1  •  ^ 

-  Pha3e{H^(ft^+ft^«-ftjl^)} 

•*’* 

Fig.  *1.  Third-order  nonlinear  system. 

This  expression  for  the  third-order  response  can  be  expressed  more 
concisely  with  the  following  notation, 

M  M  M 
^  i  s  k  S  1 
•  <=os(<»i±aktDt)t 

To  give  the  reader  a  feel  for  the  complexity  of  the  cubic  nonlinearity 
response  to  a  sum  of  sinusoids,  consider  an  input  consisting  of  the  sum 
of  three  cosine  waveforms,  with  distinct  frequencies.  In  the  real-world 
such  a  signal  might  arise  from  an  amplitude  modulated  carrier,  wherein 
the  modulation  signal  is  a  low  frequency  sinusoid;  and  might  be  the 


third-order  nonlinear  transfer  function  of  the  power  amplifier.  The 
input  is 


x(t)  =  cos(J3^  t)+A2COs(fl2t)+A^cos(n^t) 

Using  expression  (29)  the  output  is 

y(t)  -  J  ir„,|a,,,|[3|H^(i!,)|co3((ii,)t  .8  ) 

•|Hj<3a, ) |cos((3a, )t  *36^^ **3  ) ] 

.|H^(20,-a2)|<=oa((2a,-02)t 

♦r,22l«,22ll[|H<,(a,-2a2)|oo3((»,-2opt 

)  |cos((Q^  )t  +26^^  ) 

* |H/», *202) |ooa( (a, .20j)t  *2e„^.*„^ .20^) ] 

*^222 ^  ^222  ^  ^ ^^2^  ^ ^ ^2^ ^  ^ 
•|Hj(30j)|oo3((3n2)t  *380  •♦ja  1 

•r,l3l0,l3l(^l"a'“3’l°°®““3>'  *^®a  > 

.|H^(2a,.a3)|co3((2a,»a3)t  .2»„^.  »a^**2o^.a^) 
.|H^(2a,-a3)loo3((2a,-a3)i  •e2a^-Y*2a^_a^)l 
*  23 1  23  I  U  “r^2"‘’3^  I  ^  "“2“”3^  ^  *  ®a^  ~®a2~*a3*'^a^  ^ 

•l"d<“i‘V“3’l°“““i*”2'“3’'^  •9a/9a^-«„3*»a^.a2-a3> 
*l»d(‘>,-V“3’l°°’““>-V»3’‘  ‘9a,'\*®a3**a,-aj*a3> 
*l«d<»,*V‘’3’l°°9‘<”i*V”3>'^**a,*\**a3**a,.aj.a3>l 

• r, 33 1 0, 33 1 (H,( a, -283) I oo3( ( a, -2a  ) t  . e  -26a  .*a  .2a  ) 


+2|H^(Q^ )  |cos((Q^  )t  +20^^  ) 

*''223l‘’223ll^l"a‘“3’l'°*‘'“3>'  *^V' 

•|H^(2aj-a3)|cos((2a2-a3)t 


*lH^(2a2-a3)|ooa((2aj.a3)t 


•  '■233 1 "233 1  f  ‘'2*2“3’  I ‘  V2“3’ ‘  > 

+2|H^(n2)loo3((«2)t  +2ejj  ) 

*r333|Q333|[3|Ha(a3)|cos((n,)t  ) 


3'“  ''a. 


>lH^(3a3)|co3((3a3)t  +39^,  ■^♦3^,  )]}  (30) 

This  output  expression  is  only  for  the  sum  of  three  cosine  inputs;  if 
the  input  were  the  sum  of  seven  or  eight  cosine  waveforms,  one  can  see 
the  difficulty  in  calculating  the  response  in  terms  of  the  frequencies 
generated. 


Conversion  of  The  Third-Order  Volterra  System  Response 
Into  The  Third-Order  Nonlinear  Response 


The  response  of  a  third-order  Volterra  operator  can  be  expressed  in 
the  functional  form 


y^(t)  -  H,[x(t)] 


in  which  h^Ci^.T^.t^)  is  the  symmetric  third-order  Volterra  kernel  of 

the  third-order  Volterra  operator  H^.  In  this  section,  our  objective  is 

to  express  the  Fourier  transform  of  the  response  Y^(to)  in  terms  of  the 

Fourier  transform  of  the  input  XCoi)  and  the  third-order  NLTF 

.(i)^  ,10^) .  Note  that  the  arguments  are  taken  to  be  u,  rather  than  joj 

to  simplify  the  notation.  Earlier,  we  saw  how  the  transformation  of  the 
second-order  Volterra  associated  response  into  the  actual  response 
accomplished.  We  follow  the  same  procedure  for  the  case  of  the  third- 
order  Volterra  operator.  The  third-order  associated  response  is  defined 
as 

00 

y^^^(t^  h^(T^  ,T2,T^)x(t^-T^  )x(t2-T2)x(t2-T^) 

•  dT^dx^dx^  (32) 

The  third-order  response  y^it)  is  readily  obtained  by  setting  each  of 

the  arguments  of  y^^^  equal  to  t.  That  is 

y^(t)  =  y^^j(t,t,t)  (33) 

Let  us  next  compute  the  three  dimensional  Fourier  transform, 

(u)^  ,(j)2,w^)  of  the  associated  functici  y^  ^^  ( t^  ,  t^ ,  t^ ) .  Then  compute 

Y^(u))  from  this  Fourier  transform  of  the  associated  response. 

The  three  dimensional  Fourier  transform  of  the  associated  response 


is. 


^(3)^“l 
where  X(u)  i 


=•  H^(u)^  ,h)2,'D2)X(u)^  )X(w2)X(u)^) 
s  the  Fourier  transform  of  the  input  signal. 


(3M 

Taking  the 


three-dimensional  inverse  transform  yields 


r-;  ▼  .  «—  B^n 


"i  ;.”*  ’  “^  .  ^ '.  '*7  "7''.' 


7'"r  ".'  \*  ■  I  "irwi’ 


y(3)(t^  .t^.t^) 


(2it) 


L_  rrr  V  r 


T':^77T7* 

m 

.;■•  nfj*  I 


•  do)^  doj^du^ 

Setting  t^=t2=t3=t  in  our  previous  equation  we  obtain. 


y3(t) 


(.2-n) 


^  “  j(a)^ +t02+(ij3)t 

— 3  ///  3j  (o)^  ,0)2 . W3 )e  do)^do)2do)3 


(35) 


(36) 


The  inverse  Fourier  transform  of  Y3(o))  is  by  definition 


y3(t)  = 


1 


;  Y3(o))e'^“’^  duj 


(37) 


(2-n) 

Thus  we  wish  to  make  a  change  of  variable  in  the  Fourier  transform  of 
the  associated  response  that  will  give  us  the  Fourier  transform  of  the 


actual  response.  Let  o)^+o)2+o)3  =  o),  this  enables  us  to  obtain  Y^iui) 

from  Y^3^  (o)^  ,o)2.o)3).  We  shall  make  the  -desired  change  of  variable  in 
two  steps.  First,  in  eq.(38)  we  let  to  obtain 


y3(t)  = 


,  "  j(u)^+'J'^)t 

— 3  Iff  Y^3^  (o)^  ,i()^ -o)3,o)3)e  do)^  dt(/^  do)^ 


(2Tr)^  -= 


(38) 


.•  'J 


4 


.  %  %  \ 

'•*  Sv.’ 


. 


For  the  second  change  of  variable,  let  ^^+u),j=(o  in  (ilO)  to  have 


y3(t)  =  — ^ — 3  ///  Y^  ^  (u)-i|)^  ,1);^ -u)_,u)3)e'^‘^^d\j)^  dw  du  (39) 

(2Tr)^ 


Note  that  this  triple  integral  is  in  the  form  of  an  inverse  Fourier 
transform  of  dimension  one.  From  this  observation  we  can  rewrite  it  as 


' ' '  ^ 


Y^Cu) 


1 


(2it)  - 


2  7  /  Y^3j  (a)-i((^  ,  -u)3,u)3)dii)^  do)3 


(2tt) 


2  /  /  Y^3^  (u)-4i^  ,\|)^ i|)2)dij)^  dij)2 


(^0) 


(41  ) 


v 

vvSS 
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In  other  words,  Y^(io)  is  the  integral  of  3)  ('('i » '('2’ ’*'3^  along  the  plane 

w  i  alternatively  Y,(oj)  is  the  profile  of  Y,_.  (i|;,  along 

I  j  i  ^5)1^5 

the  line  =(1)2=41^. 

CASCADING  OF  NONLINEAR  SYSTEMS 

Before  we  deal  with  the  matter  of  compensation  of  nonlinear 
systems,  we  must  first  characterize  a  tandem  connection  of  two  general 
nonlinear  systems.  Fig.  5  shows  the  cascade  connection  of  such 
Volterra  systems.  In  this  section  we  wish  to  derive  an  equivalent 
Volterra  NLTFs  description  of  the  two  nonlinear  systems  in  tandem. 
Suppose  that  the  two  subsystems  are  characterized  by  the  NLTFs  and 

P^,  n  =1 ,2  , . . . ,  ®  . 

Then 

00  00 

w  =  H[x]  =  I  H  [x]  =  I  w  (42) 

n=1  n=1 

and  the  output  y  is  given  by 

00 

y  =  p[w]  =  I  P  [w]  (43) 

k=l 

00  00 

-  I  I 

k=1  n=1 

To  use  the  nonlinear  current  method  [3],  let  us  replace  x  with  ex  to 
form  a  new  output; 

00  00 

y  -  I  Pj  I  H^Cex]].  (44) 

k-1  n=1 

■  I  pj  i  .\cx]i 

k-1  n-1 
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•:%w 


> 
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Fig.  5.  Cascade  of  two  nonlinear  systems  H  and  P. 

Now,  from  the  basic  k-linear  property  [1]  of  the  operator  Pj^,  it  is  seen 
that 

n.  n. 


Pj.[  I  ”  I  •••  I  e  ^••*g  ^  P^^^n  ’’“’''n  ^ 

n-1  n^-1  n^^-l  1  k 


('♦S) 


in  which  P,^l'}  is  the  k-linear  Vol terra  operator 

''kl'n,  •••••".  I 

1  k 


00  00 


///•••/  '**  ’^n  *•* 

1  k 


—  00  —00 


•  -  w  *  » 


L*  *  •  ^ 

V  V  4 
ym  •■ 


.•.“-.■.'I.' 

rur. 


Then 


m  <”  .., 

y-l  Z  •••  I  c  ' 

m-1  n  -1  n_-1 

I  m 


09  n.  •'■•••+n 


P  |w  ,  •  •  • ,  w  I 

m'  n  ’  ’  n  ' 

1  m 


{>i7) 


This  relationship  can  be  re-expressed  in  terms  of  the  ”x-y  transfer 
function"  Q  as 
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9.*^  •' 

•.%C- 


Cs.-. . . • . .  • . . .  . 

J^’v  v  vS  '  -  V  -'  ■ 


-'In'. 


^  •  N  *  te 

•V  .V' 


00 


y  =  I  QM  (‘»8) 

k=1 

Now  if  like  powers  of  e  are  collected,  the  Q’s  can  be  determined  using 
"operator  notation"  as 


s  = 

PlH, 

«2  = 

+P  H 

2  1 

s- 

P1H3 

>2P2{H^.H2} 

+P  H 

31 

PlH, 

+2P2{H^.H^} 

+P2CH2]  +3P3{H^.H^,H2)  +P^[H^] 

S  = 

^«5 

>2P2{H^.Hj^) 

+3P3{H^ 

Q  =  P.H  +  terms  involving  P.  for  j  S  (m-1)  (50) 

m  1  m  j 

This  relationship  is  needed  in  our  discussion  of  the  mth-order  nonlinear 
inverse. 

These  equations  can  be  expressed  in  the  Laplace-domain  as 


Q^(3)  »  P^(s)H^(s) 


(s^ +S2)H2(s^  ,S2)+H^  (3^  )H^  (S2)P2(S^  .s^) 


P^  (s^+S2*3^)H2(s^  .s^.s^)  +  3  [H^  (3^  )H2(S2,3^)P2(s^  ,S2+3^) 
(s2)H2(3^  ,s^)P2(32,3^+32)+H^ (s^)H2(s^ , s^ )P2( ,3^ +3^ ) ] 


(s^  )H^ (s2)H^ (s2)P2(s^ ,3^.3^) 


(51  ) 

Our  objective  is  to  reduce  the  nonlinear  response  in  a  system  without 
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I 


distorting  the  linear  characteristics.  Since  the  linear  response  must 


not  be  changed,  the  form  of  should  be  ;  and  therefore, 
^ 


(52) 


where  I  is  the  identity  operator.  Furthermore,  the  second  and  third- 
order  net  nonlinear  transfer  functions  (Q^  s^d  Q^)can  be  made  zero  by 
choosing  and  as 


P2  =  -P^H^H^ 


-1 


(53) 


P^=  -(H^  ^2P2{H^,H2})H^ 


-1 


Example 

A  simple  example  is  used  to  show  how  a  nonlinear  compensator  is 
constructed.  Only  the  second-order  distortion  is  reduced  in  this 
example,  but  in  Chapter  4  we  will  address  both  second  and  third-order 
nonlinearly  generated  distortion.  Suppose  the  NLTFs  of  a  circuit  with  a 
single  nonlinear  device  [3]  are: 


(s) 


Cs 


Ls 


Cs  Z(s) 


(54) 


LCs  +Lg^s  +1 


H2(s^,32)  =  -a^Z(s, +s^)  {H, (s, )H, (s,)} 


1  ''Z' 


1 '1  '1 '2' 


Then  becomes 


^2^®1’®2^  =  -H2(s^ (s^ +3^) 


=  a 


LC(3^  ^s^) 


^  LC(s^+S2)^  +Lg^  (3^+82)  +1 


(55) 


The  algebraic  expression  for  P^(s^,S2,s^)  is  determined  by  evaluating 
the  second  expression  of  equation  (55).  While  this  evaluation  is 
mathematically  straight  forward,  the  resulting  expression  is  involved. 
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Therefore,  the  specific  form  of  will  not  be  determined.  This  then 
completes  this  example. 

THE  APPLICATION  OF  VOLTERRA  THEORY  TO  REDUCE  NONLINEAR  EFFECTS 

Volterra  theory  may  be  used  to  predict  and  diminish  undesirable 
effects  of  the  nonlinearity  of  a  system.  One  possible  approach 
suggested  by  Wiener  [1]  is  to  use  linear  feedback  to  reduce  the 
nonlinear  contributions.  Let  us  review  the  feedback  method  briefly. 

The  reduction  of  nonlinear  response  by  linear  feedback  was  applied 
to  transistor  amplifiers  by  Narayanan  [5],  In  that  paper  the  open-loop 
transfer  function  of  a  transistor  amplifier  is  characterized  by  a  sum  of 
Volterra  kernels,  shown  in  Fig.  6.  These  kernels  take  into  account  the 

4 

frequency-dependent  nature  of  the  nonlinearities. 

The  basic  closed  loop  equations  for  the  circuit  are: 
e  -  x-v 
y  =  H[e] 

V  =  B[y] 

y  =  GCx]  (56) 

The  closed  loop  system  operator  G,  expressed  so  that  it  satisfies  all 
input  signals, 

G  =  H  [I-BG]  (57) 

where  I  is  the  identity  operator.  The  closed  loop  kernels  for  the 
first,  G^(a)),  second  ,0)2 ) .  and  third-order,  G^ ( u)^  , 0)2 , o)^ )  are 

obtained  by  successively  equating  the  linear,  second  and  third-order 


terms. 


pryjyyrrv  r^.r-  r--..  -..y  v  » V  V  >.■  4»IT'^  *;.*■.'  rj  ■'»  WJ  WIW.'WVW. 


I 


c 

V 


^2^“l ’“2^ 


H2(o)^  1(1)2  ) 


[  1  ((0^  )B(  (D^  )  ][  1  +H^  ((i)2)B((D2)  ]  1  (o)^  +0)2  )B((d^  +0)2) 


(59) 


G-((d^  ,(1)21(1)2)  -  [H2((i)^  ,(1)21(1)0) 


i»1  1 +H^  ((Dj^)B((D^) 


{3H2((i)^  ,(i)2+(Do) 


(u>,  )B((i.^  ) 


1b(  oj^  +(1)2  )H2(  (1)21  (112)^2  (<*>21(1)2) 


1+H^((D2)B(ui2) 


)  •  (- 


1  +H^  ((>>^)B(u^)  1  +H^  ((J2'*’(D2)B((D2''‘(D2) 


1  +Hj  (d)^  +(i)2+(D2  )B(  u)^  +(i)2'*'id2  ) 


-)] 


(60) 


Where  is  unsymmetrized.  The  application  of  feedback  has  reduced  the 

gain  of  the  linear  transfer  function  by  the  factor  L((o)»  1 +H^  (a))B(u)) . 

The  second-order  expression  reveals  that  the  application  of  feedback 
reduces  the  gain  by  the  factor  of  L()d^  )L(  0)2  )L(  (1)^ +(*>2)  •  The  third-order 

term  indicates  that  the  feedback  has  reduced  its  gain  by 

L((i)^  )L((D2)L.((D2)L{(d^ +(D2+(i)2)  (and  further  terms  at  the  frequencies  of  the 

two  possible  second  order  kernels). 

From  the  above  discussion,  several  observations  can  be  made. 

First,  the  attenuation  in  each  output  term  is  related  to  its  generating 

kernel-order;  The  higher  order  terras  are  reduced  significantly,  however, 

the  reduction  achieved  for  the  lower  order  terms  may  not  be 

sufficient  for  certain  applications.  This  approach  is  only  useful  to 

systems  in  which  feedback  can  be  provided  conveniently.  Clearly, 

-  23  - 


V.  .N  .%.%1 


iiL^ 


'-t 


.---1 


k  *  .  «  .  t  .  • 


w'-  A 


"JT' 


-V-'. 


-• 


.•  :i 

■  /■'■y 


G 


Fig.  6.  Volterra  model  of  a  transistor  with  linear  feedback. 

feedback  is  not  easily  Implemented  if  the  nonlinearity  of  the  entire 
link  —  transmitter,  channel  and  receiver  —  is  to  be  compensated. 
Finally,  using  linear  feedback  to  reduce  interference  from  the  NLTFs 
also  reduces  the  gain  of  the  linear  TF. 

In  our  method  of  post-compensation  (to  be  discussed  later)  these 
disadvantages  are  alleviated  to  a  large  degree. 
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CHAPTER  III 


INTERMODULATION  AND  CROSS-MODULATION  IN  NONLINEAR  CHANNELS 

AND  AMPLIFIERS 

The  phenomena  of  intermodulation  and  cross-modulation  occur  when 
two  or  more  signals  are  mixed  in  a  nonlinear  element  in  such  a  way  that 
multiplicative  combinations  of  these  signals  are  produced.  New 
frequencies  are  generated  (Intermodulation),  and  intelligence  from  one 
carrier  may  jump  on  to  another  carrier  (cross-modulation).  This  occurs 
in  receiver  input  circuits,  as  well  as  in  transmitters,  and  sometimes  in 
a  nonlinear  element  in  the  channel.  As  far  as  transmitters  are 
concerned,  the  process  involves  the  reception  of  an  unwanted  signal  by 
the  transmitting  antenna,  which  conducts  it  back  to  the  final  stage  of 
the  power  amplifier,  where  it  is  mixed  with  the  transmitted  signal.  The 
process  is  therefore  of  greatest  significance  when  both  the  unwanted 
signal  and  the  nonlinear  product  are  within  the  passband  of  the  final 
amplifier.  This  can  only  occur  if  the  nonlinear  device  has  a 
characteristic  of  odd  order.  In  this  chapter  several  types  of  nonlinear 
distortion  will  be  discussed. 

Distortion  Caused  by  Static  Nonllnearltles 

Static  nonlinearities  may  be  characterized  by  a  power  series  with 
constant  coefficients.  For  example,  a  single  input,  single  output  static 
nonlinearity  may  be  characterized  as 


(61) 


y(t)  =  a^x(t)+a2X  (t)+a^x^(t)+  ...  +a^x  (t)+  ... 

» 

=  J  a  x'^(t). 
n=1 

This  type  of  nonlinearity  is  sometimes  referred  to  as  a  memoryless 
nonlinearity,  which  is  to  say  that  the  present  output  is  only  a  function 
of  the  present  input  values  and  does  not  depend  upon  the  past  input 
values.  A  nonlinear  resistor  is  an  example  of  a  static  nonlinearity. 

CROSS-MODULATION  DISTORTION 

Cross-modulation  occurs  when  the  modulation  information  of  one 
carrier  migrates  to  another  carrier  that  did  not  have  these  sidebands. 

To  illustrate  the  phenomenon  that  takes  place  in  cross-modulation, 
consider  a  nonlinear  device  that  has  an  input-output  characteristic 
specified  by  a  power  series.  Let  the  input  to  this  device  be  the 
sum  of  an  unmodulated  carrier  with  frequency  and  a  carrier  with 

frequency  ,  modulated  with  a  tone  of  frequency 

x(t)  =  V( 1 +b  cos(n  t))cos(Q, t)  +  V  cos(n.t)  (62) 

1  m  1  2 

The  spectrum  of  this  signal  is  shown  in  Fig.  7.  Let  us  assume  the  terms 
of  the  static  nonlinearity  have  negligible  contribution  beyond  the  cubic 
term.  The  terms  of  interest  lie  at  the  frequencies  £J^  and  these  are 

the  cross-modulated  terms.  The  linear  part  of  the  response  is 

y^(t)  =  a^  v[( 1 +b^ cos(n^t))cos(n^ t)+cos(n2t)  ].  (63) 

Using  the  identity 

co3(a)co3(6)  =  [cos(oi-B)+cos(a+B)  ] 

(63)  may  be  written  as 


y^(t)  -  a^v[cos(n^  t)+  —  cos(n^-Q^)t  +  cos(n^ +n^)t+co3(n2t) ]  (S'*) 


G  2  CJ  G  “■ 

Fig.  7.  Spectrum  of  quadratic  nonlinearity 
generated  cross-modulation. 

The  output  spectrum  is  similar  to  (in  fact  proportional)  to  the  input 
spectrum.  The  quadratic,  or  second-order,  response  is 

y  (t)  =  J  {a  V^[2(1+cos(2n,t))+2b, (cos(n  t)+cos(2n  -n  )t) 

^  H  £1  lira  Ira 

+  2b^ (cos(n^t)+cos(2n^±n^)t)+2(cos(n^+n2)t) 

^  2 
-r-  (iicosCa  -n  )t)+b, {i+cos(2n -2n  )t) 

^  >  d  \  I  m 

+  b^(cos(2fi  t)+cos(2fi,  t)+b, (cos(n„-B  +n  )t 
I  m  11  2  I  m 

+  2co3(n  +a  -a  )t)+  -g-ci +cos(2a  +2a  )t) 
dim  d  1  m 

♦  2b^ (C03(a2+a^±a^)t*(1+  cos(2a2t))]}  (65) 

Let  us  introduce  a  notation  to  reduce  the  length  of  the  expressions. 
Let  cos(a)t  be  represented  by  and  co3(a,  +  a-,)t  by  C„  „ 


.  Then 


the  previous  expression  becomes 


y^Ct)  -  1  C 

m  m  3  ID  3 


+  2(C  +  C  )+b  (C  +C  ) 

^  a-Q^  r  2n -n„  2a +n„'' 

3m4tn  imim 


'1 


+  2(1+  -:r^)C^  +2b  (C  +C  ) 

^  2  ^  20,  K  n,-a  Q,+n 

1  4  m  4  m 


^^i^^2a  -2n  *Sa  +2a  ^*^2a 

Imim  2 


(66) 


even  more  compactly, 


y,M  .  J  (a,v^[2(2,  .afc,„  •'•c  .2b,c 

m  m3 


3"  m 


*  — ’'2a, .a -"Cb  •"'=20 

1  1  m  M  4  m  2 


(67) 


where  »  ^2~^1  ^4“  ^2*^1 


Notice  that  there  are  no  cross-modulation  terms  in  this  case. 

The  cubic  or  third-order  terra  in  the  series  requires  a  method  for 
tabulating  the  frequency  components.  This  will  allow  us  to  view  the 
various  terms  faster  and  this  approach  can  be  used  to  analyze  higher 
order  nonlinearities.  For  the  input 

b. 


'1 


x(t)  =  V{C  +  — ^(C  +C  )  +  C  } 

'  ^  ^  fl,  2  a  -a  a  +n  a. 

1  1  m  1  m  2 


the  cubic  output  is 


y^(t)  -  a^x^(t). 


We  may  write  the  input  in  an  expanded  form,  namely 
x(t)  =  A.x,  (t)+A-X_(t)+A  x_(t)  +  A,  x,,(t) 


=  A^cosCSj^t) 

Then  the  output,  y^Ct)  may  be  computed  by  a  systematic  combination  of 

the  Input  constituents  as  follows:  the  order  is  three  therefore  the 
number  of  columns  (containing  the  input  constituents)  is  also  three. 

See  Table  1.  Each  row  represents  a  distinct  combination  of  the  three 
constituents  from  the  available  four  sinusoids.  For  y^(t)  the 

components  of  the  output  are  listed  in  Table  1.  The  multipliers  1,3  or 
6  in  the  right  half  of  the  table  arise  from  the  number  of  distinct 
permutations  of  the  Indices  of  the  constituent  sinusoids. 


Using  the  cosine  identities 


003(a)C03(6)003(»)  .  5 


cos  (a)Cos(e) 


aos^ia)  -  ;  (3C  •C,  ) 

4  a  3a 


The  terras  that  produce  cross-modulation  are 


3V^b.C  C  C  +3V^b.C  C 

1  1  m  2  1  m  2 


?  b  V^{C  +C  + 
H  n -R_  n 

ra  2  m2 


♦C  +C  } 

-R  +R_  -R  -R, 
ra  2  m2 


=  ^  b,V^{C  +C  )■*■ 
2  1  R--R„ 

2  m  2  m 


Actually  since 


"3  ” 


3a,b, V^co3(R_t)co3(Q„t)  +  ... 


TABLE  1 

COMPUTATION  OF  OUTPUT  FREQUENCIES  GENERATED  BY 
THIRD-ORDER  NONLINEARITIES 


permutations 
1  1  1 

1  1  2 

1  1  3 

1  1  H 

1  2  2 
1  2  3 
1  2  4 

1  3  3 

1  3  4 

1  4  4 

2  2  2 

I 

2  2  3 

2  2  4 

I 

2  3  3 

2  4  4 

i 


translation 
A^cos^Cn^  t) 

3A^A2  cos^  ( t )  cos  ( £22 1 ) 

3A^A2Cos^(£2^  t)co3(n2t) 
2  2 

3A^A^cos  t)cos(nj^t) 

3A^A2C03(f2^  t)cos^(n2t) 
6Ai  A2A2C03(£2^  t)cos(£22t)cos(n^t) 
6A^  A2A^C03(12^  t)cos(£22t)co3(nj^t) 

SA^AjCOsCD^  t)c03^(£22t) 

6A^  A^A^^co3(£2^  t)cos(£2^t)cos(B^^t) 

2  2 

3A^  A^cos(n^  t)cos  (Rjjt) 

A2C03^(£22t) 

3A2A^co3^(a2^)co3^^3'^) 
2  2 

3A2A^co3  (£22t)co3(nj^t) 

3A2A^co3(R2t)co3^(n^t) 

2 

3A2Aj^co3(n2t)cos  (flj^t) 


V  Vv 


S.-  K.--J 


-v 

.v-v^^yo 


-.V.V. 


•  •  •.  •  •  ^ 


/■■•'.■•-.'I 


■  ■  1 
■  ■  •i 


v.'.- -J 


3  3  3 


30  - 


A^cos^Cn^t) 


Table  1  (cont'd)" 

2  2 

SAjA^cos  (U^t)oo3(:i^t) 

.  W'-  ' 

3A^A^co3(f22t)cos^(£lj^t) 

A^cos^(n„t) 

The  carrier  term  may  be  approximated  as 


a^V  003(^2^)+  ^  a^V^co3(n2t) 


a^  VcoaCn^t) 


Thus,  the  cross-modulation  at  <^2  is 


y^„(t)  =  a,V(l +b^„co3(n^t))co3(«,t) 
cm  1  cm  m  2 


.V_‘. 


3  2 

b  =  3  -r^  b.v'" 
cm  ^ 


The  cross-modulation  factor  is  the  ratio  of  b  to  b,  is 

cm  1 


cm 

cm  =  — ^ 
o  b, 


As  an  example  consider  a^  =  10,  a^  =  0.02,  b^  =  0.5,  V  =  1  volt.  Note 

that  the  amount  of  modulation  is  50? .  Then 

b  „  =  3(0.002)(0.5)(1)^ 
cm 

=  0.003 


V  (t)  -  10(1+  0.003co3(n  t))co3(n_t) 
cm  m2 
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A'.‘  •  ■ 

.''"I 


0,006 


^  tr-  V'-  V'-  «ri  ■ 


CM 


cm 


'1 


This  index  measures  the  amount  of  modulation  which  has  leaked  onto  the 
unmodulated  carrier. 

CROSS-MODULATION  FOR  NONLINEARITIES  WITH  MEMORY 

Let  the  input  be  the  sum  of  a  continuous  carrier  and  an  amplitude 

modulated  carrier  as  before 

x(t)  -  Vd+b.cosn  t)co3n.t  +VcosD_t 
1  m  1  2 

this  can  be  written  in  a  complex  phasor  form  as 


x(t)  =  [  2e  +b^  e 


+D^  e  +2e 


~jn. t  -j(n  -n  )t  -j(n.-n^)t  -jn.t 

>2e  1  .b^  e  ^  ^b^  e  ^  m  ^^e  ^  ] 


(72) 


To  reduce  complexity,  let 


e  =  eJ^^ 


Then  we  may  write  x(t)  as 


x(t)  .  j  [  2e  e  *6,  v.„  3  •  CC 

1  1m  1  m  2 


(73) 


where  CC  stands  for  complex  conjugate  terms. 

If  we  use  the  tabulation  technique  shown  earlier  and  select  only  the 

terms  that  contribute  to  the  cross-modulation  namely,  the  output 

comLonent  at  the  frequency  n„±n  turns  out  to  be 

2  m 

2  m  2  m 

+  H^(-fl^  ’  ~”r^m’ ”2^®n2"^ni 


CC] 


(74) 


To  discern  the  elements  that  cause  cross-modulation,  we  need  to  express 
this  equation  differently.  We  will  assume  is  a  continuous  function 

and  that  the  frequency  of  modulation  is  much  less  than  the  carrier 
frequency, 

m  1 

We  may  now  write  the  cross-modulation  terms  generated  by  the  third-order 
nonlinearity  as 

+£i  ^  -n 

^  ta  2  m  2  m 


=  2  K  v3 


I  b^V^|H2(a^,-0^,fi2)|Cco3((fi2-0n,)t 


♦cos((a_+n„)t  +B  ] 
2  m  2 


where  62=  phaseCH^Ca^ ,-fl^ ,«2)) 


If  we  look  at  the  third-order  system  response  to  the  carrier  term 


^3  ’  2  ^/‘^2^®a2*  2  *^3^”2’”2’~”2^®a2"^ 


2 


=  v|h^ ( a^) I cosCn^t*  ) 


where  =  phase {H^(a2) 


Thus  the  a^t  term  becomes 


■  V|H,(a,)(cos(a2t  *6,) 


+3  b, I ( a^ , -a, , a^ ) I cos( a^  t  +  B2 ) cos ( a^t ) 


H,  (a^ ) I V[co3(a2t  +6^ )+b^cos(a2t  +B2)co3(a^t)]  (77) 


Therefore, 

-  |H  (Q  )  |v|[co3(fl  t+B,  )+b  co3(a-t+6_)cosa„t]  (78) 

cin  '1^''  ^ix  d  d  n 

The  cro33-modulation  term  I3  different  from  the  memoryle33  caee,  there 

wa3  pure  amplitude  modulation.  Expresaing  the  croaa-modulation 

content  as  a  function  of  6^  and  we  must  consider  two  cases: 

Case  I 

ycn,(t)  -  |H^  (fl^)  [Vd +b^co3Q^t)co3(n2t+B^ ) 

This  case  generates  pure  amplitude  modulation. 

Case  II 

62  =  61  -  1 

then  the  cross-modulation  is 

where  b  <<  1 

X 

This  case  generates  pure  phase  modulation. 

INTERMODULATIOM  DISTORTION 


Intermodulation  distortion  is  the  process  whereby  beat  signals  at 
various  sum  and  difference  frequencies  are  produced.  These  signals  are 
Impossible  to  filter  with  a  linear  filter  because,  the  Interference 
frequencies  lie  within  the  passband  of  the  uncorrupted  signal.  For 
example,  consider  an  Input  consisting  of  the  sum  of  two  cosine 
waveforms,  _  ^4  _ 


w*v«r.-wvirj 


»rf.'f  .  vj’T'v  V  '.n- 


x(t)  -  VC  +VC 
“l  “2 

The  intermodulation  terms  that  would  be  generated  by  a  second-order 
nonlinearity  acting  on  this  Input  are 

The  third-order  intermodulation  products  are  at  the  frequencies 
(2n^±a2)f  (202*^1)  and  could  lie  within  the  passband  of  either  carrier. 

The  complete  expression  is  found  to  be 


^3iM  “  i  ±a2^'^2a2±ai  ^ 


(79) 


■  If"!  M'  VW) 


•_  *-  •. 


The  magnitude  of  these  IM  terms  is 
D, 


2  a.v3 


'3  ~  »»  3 
Then  the  IM  index  is  expressed  by 


(80) 


D3 

IM3  .  ^  V 


3a 


^a 


3  ^ 


(81) 


1  "“1 

This  concludes  our  discussion  of  interference  types. 


CHAPTER  IV 


COMPENSATION 


Communication  channels,  both  microwave  and  satellite,  generally 
exhibit  nonlinear  characteristics.  When  frequency  division  multiplexing 
is  employed,  it  is  Important  to  carefully  assess  the  extent  of  possible 
intermodulation  distortion.  When  this  distortion  exceeds  permissible 
limits,  it  is  reasonable  to  ask  if  it  can  be  reduced  by  suitable 
compensation.  In  this  chapter  we  present  the  compensation  of  both 
second  and  third-order  NLTFs.  The  design  methodology  developed  here  has 
been  implemented  in  FORTRAN  programs  VC0MP2  and  VC0MP3  [12],  [13]. 
Results  of  example  runs  on  these  programs  are  presented,  which  confirm 
the  effectiveness  of  our  technique  in  reduction  of  IM  interference  in  a 
specified  band. 

The  Pth-Order  Inverse 

Suppose  the  input-output  relationship  of  Fig.  1,  representing  a 
weakly  nonlinear  system,  is  given  by 
y(t)  »  H[x(t)] 

-  H^[x(t)]+H2[x(t)]+H^[x(t)]+  •••  (82) 

then  a  p-th  order  post-inverse  is  defined  as  follows.  A  nonlinear 
system  G  (also  assumed  to  be  representable  in  a  Volterra  series)  is  the 
p-th  order  post-inverse  of  the  nonlinear  system  H  if  the  composite 
system  operator  Q,  given  by. 


Q[x(t)]  -  y(t) 

-  G[H[x(t)]]  (83) 

contains  only  nonlinear  terms  of  order  p+1  or  higher  and  Its  linear  term 
Is  the  Identity  mapping.  That  Is 

(84) 


QCx(t)]  -  x(t)  +  Q_,[x(t)]  +  ••• 

p+1 


In  this  report  our  concern  Is  limited  to  only  second  and  third- 
order  nonllnearltles.  We  give  below  the  conditions  upon  the  NLTF's  of  H 
for  the  cases  of  second  and  third-order  post-inverses:  The  following  two 
conditions  must  be  satisfied  [1]  for  a  second-order  inverse  to  exist: 


(85) 

'1»2^^2«1  =  0 

(86) 

For  G  to  be  the  third-order  inverse  of  H,  three  conditions  must  be 
satisfied,  namely  (85),  (86)  and  the  following 
G^H^  -  G2{H^,H2}  *  G^H^ 

(87) 


Note  that  if  »  0,  l.e.,  the  second-order  nonlinearity  is  absent,  then 


(87)  reduces  to 


G^H3  >G3H,  =  0 


Although  the  concept  of  a  p-th  order  inverse  is  a  useful  one  from  a 
theoretical  standpoint,  its  practical  design  utility  is  quite  limited. 

In  fact,  even  the  first-order  inverse  defined  formally  by  (85),  has  a 
somewhat  limited  utility  in  practice.  Indeed,  such  an  inverse  would 
require  the  following  relationship  to  be  satisfied 


*  -  • 


•  * 


mm 


m- 


;.\v 


G,(s)H,(s)  =  1 


(88) 


Implying  that  G^(s)  be  the  inverse  of  H^(3).  This  is  impractical 

for  two  reasons.  First,  this  would  necessitate  infinite  gain  for  G^(s) 

at  values  of  s  where  H^(s)  is  zero;  second,  if  the  transfer  function 

H^(s)  has  zeros  in  the  right-half  plane,  then  G^is)  would  contain  poles 

in  the  right-half  plane  and  would  thus  be  unstable.  In  most 
communications  applications,  a  global  inverse  is,  in  fact,  not  required. 
Indeed,  if  the  signal  spectrum  is  limited  to  the  band  then  the 

objective  should  be  to  design  G^  so  that  G^H^  may  appear  to  be  almost 

equal  to  one  (or,  more  generally,  equal  to  e  in  this  band. 

As  stated  above,  application  to  communication  systems  requires 

minimal  transmission  distortion  only  over  a  particular  frequency  band 

(or  bands).  Let  us  denote  this  band  (or  union  of  bands)  as  <R  the 

roi 

subscript  'roi'  signifies  the  region  of  interest.  Also,  let  us  denote 
the  desired  characteristic  of  the  composite  system  as  C(u).  Then,  it  is 
reasonable  to  state  the  desirable  attributes  of  the  compensator  as 
follows;  1)  The  composite  system  Q,  should  have  a  characteristic  as 
close  to  the  desired  as  possible,  i.e.,  should  approximate  C{u) 

over  the  band  <R  ,>,  2)  the  contributions  of  the  NLTFs  H_,  etc.,  in 

the  band  ^nd  possibly  other  interference  regions  of  interest 

should  be  minimized,  3)  the  compensator  should  be  stable,  and  ^) 

the  structure  of  the  compensator  should  be  simple  and  relatively 
robust,  i.e.,  insensitive  to  minor  changes  in  the  parameters  of  ,  H^, 

etc. 

In  the  Appendix  we  present  two  simple  examples  dealing  with  the 


compensation  of  linear  distortion  in  bandpass  systems.  These  examples 
are  useful  in  that  they  provide  Insight  into  the  structure  of  the 
compensator  TF, (compensator  pole  location  with  respect  to  the  system 
bandpass  characteristic)  needed  to  reduce  linear  distortion.  This 
provides  the  basic  compensator  TF  form  to  reduce  the  distortion  caused 
by  NLTFs  in  band-pass  systems. 

COMPENSATION  OF  THE  (KIADRATIC  SYSTEM 

A  communication  channel  that  exhibits  only  quadratic  nonlinearity 
(and  no  higher  order  nonlinearity)  may  be  represented  by  a  linear 
transfer  function  in  parallel  with  the  quadratic  TF,  as  shown  in  Fig.  8 
This  model  might  represent  a  transmitter,  the  channel,  the  front-end  of 
a  receiver,  or  a  combination  of  these  three. 

When  an  input  signal  consisting  of  the  sum  of  two  sinusoids, 
x(t)  -  V^cosfl^t  +7200302^ 

is  applied  to  the  system  of  Fig.  8,  the  output  is  the  sum  of  the  linear 
and  quadratic  responses.  The  response  of  the  linear  block  is 
w^(t)  -  (n^ ) lco3((a^ )t  +0jj) 

>V2|H^(n2)|c03((n2)t  ♦Ojj)  (89) 

where 

=  PhaseCH^  (fl^  ) ) 

The  response  of  the  quadratic  block  is 
W2(t)  = 

^  [v^|H2(£3^,n,  )|cos((2n^)t  +^,)+V2lH2(n2.«2H°os((2fl2t 
+2V^V2|H2(n^ ,-n2) 1co3((Qt-«2)^  **3^ 


+2V^  ,fi2)  +82 )t 


+2V^H2(R^ ,-n, )>2V2H2(Q2.-fi2^ ] 

where 

(|)^  -  Pha3e{H2(n^  ) } 

(ti^  =  Pha3e{H2(fl2*^2^  ^ 
ij)^  =  Pha3e{H2(n^  ,-n2) ) 

())j^  -  Pha3e{H2(n^  ,^2^} 


(90) 


LINEAR 

OPERATOR 


QUADRATIC 

OPERATOR 


[> 


Fig.  8.  Nonlinear  communications  channel  model. 


The  system  output  then  contains  components  with  the  following 
frequencies:  ,  02'  20^,  2021  S  ”^2 

For  reasons  that  will  shortly  become  apparent,  consider  that  0^  and 
^2  ai'e  nearly  equal.  Thus  as  shown  in  Fig.  9  the  intermodulation  terms 


occur  near  the  baseband  and  around  the  frequency  20- . 


In  realistic 


situations  the  input  spectrum  is  generally  distributed  over  a  frequency 


band  with  a  bandwidth  of  r  rad./s.  This  type  of  input  might  represent  a 
communications  signature.  The  regions  of  the  output  spectrum  where  IM 
distortion  occurs  may  be  obtained  by  substituting  the  boundary  values  of 
the  input  frequency  band  for  and  from  the  previous  discussion, 

“l  ■  2  •  “2  ■  “c-  2 

it  is  readily  seen  that  the  output  spectrum  would  now  be  as  shown  in 
Fig.  10. 


Fig.  9.  Spectra  of  the  input  and  output  of  the 
quadratic  nonlinear  system. 

Note  that  the  output  interference  energy  is  spread  over  two  separate 
bands  (one  near  the  zero  frequency  and  the  other  at  double  frequency). 

A  communications  system  may  be  compensated  by  pre-distorting  the 
input  signals  with  a  nonlinearity  which  is  an  approximate  'inverse'  of 


the  channel  nonlinearities.  Alternatively,  one  may  use  post¬ 
compensation  wherein  the  compensator  is  placed  after  all  of  the  system 
nonlinearities  as  shown  in  Fig.  11,  for  example,  after  the  front-end  of 
a  receiver.  Bell  Laboratories  have  used  pre-distortion  for  compensation 
on  a  narrow-band  (compared  to  the  carrier  frequency)  signal  [6].  Their 
compensator  structure  as  well  as  the  theory  developed  are,  however,  not 
adequate  for  more  general  applications.  We  will  use  a  post-compensation 


INPUT 


Fig.  10.  Spectrum  of  the  multiple  frequency  input  to  the 
quadratic  nonlinear  system. 

strategy  and  will  apply  Volterra  theory  to  develop  a  general  approach 
for  compensator  design.  With  minor  modifications  the  methodology  can 
also  be  extended  to  pre-compensator  design. 

Post-Compensator  Representation 

The  post-compensator  transfer  characteristic  should  have  sufficient 
latitude  so  that  by  proper  parameter  selection  reduction  of  the 

-  kl  - 


5 


nonlinear  effects  is  possible.  The  post-compensator  structure  shown  in 
Fig.  12  appeared  promising  and  was  proven  effective  in  our  studies.  The 
sections  of  this  compensator,  are  a  gain  constant  for  the  linear  kernel 
and  a  parametric  nonlinear  block.  Only  the  nonlinear  response  of  the 
System  Under  Compensation  (SUC)  will  be  compensated;  the  linear  response 
will  be  left  essentially  unaltered.  The  output  frequencies  generated  by 
the  channel  NLTF  are  out-of-band  Interference  terras.  Our  compensator 
will  reduce  the  interference  at  the  frequency  2n^ .  Having  chosen  the 

interference  band  of  interest,  the  compensator  design  will  now  be 
discussed.  The  output  of  the  compensator  is 

y  =  yi-y2 

where  the  output  of  the  linear  TF  of  the  compensator  is 

y^  =  [w^ ]+G^ [w^]  +  •••  (91) 

The  compensators  quadratic  output  is 

where  is  the  bilinear  Volterra  operator  [1]. 

Since  our  analysis  is  restricted,  for  the  moment,  to  the  second-order 
response,  we  will  omit  all  terms  shown  underscored  in  (91)  and  (92). 

Now,  substituting  for  w^  and  w^  from  (89)  and  (90)  and  recognizing 

that  G^ (s)  =  K, 

y^  =  K[vjH^(q^)|cos((n^)t  +ejj)+V2lH^(n2)|cos((n2)t 

+  ^  [V^|H2(n^  )|cos((2n^  )t  ^^2^ 

.-p^)  |cos((n^-n2)t  +(J.,,)+2V^V2|H2(n^  ,^2)  I 

•  cos((n^+n2)^  +4>j,)+2V^H2(ni  ,-n,  )+2V2H2(n2-"fi2^^^ 


(93) 


SYSTEM  UNDER  COMf»ENSATlON 


COMPENSATOR 


where  K  is  a  scalar  gain.  The  quadratic  response  is  given  by 

^2  °  2  tvJlH^(n^)|^lG2(Q,.Q^)|cos((2n^t)+2ejj^+<D^) 


*2[v^v^rlH^(c}^)H^((2pc^(a^,-(iplco3((i}^-(}^)t  +  *4,^) 


*2[v^v^rjH^(si^)H^(!i^)G^(a^,a^)jco3(({i,*[i^)t  *  e„  +e„  +(^„) 


1  “2^"  “^2 


+2vjH^(n^  )G2(fi^  .-n^  )+2V2H^  (f22)G2(n2,-R2)] 

Note  that  the  right-hand  side  of  (9M)  represents  G^Cw^ ]  since  all  other 

terms  in  (92)  were  omitted.  From  the  earlier  section  we  know  that 
removal  of  the  second-order  output  effects  requires 
GiH,  =  I 

S«2*°2»1  “  ° 

or,  in  terms  of  the  composite  system  TFs 

S  =  I 

Q2  =  0 

Since  we  stated  that  our  goal  is  to  reduce  the  distortion  at  2R^ ,  we  do 
so  by  setting  terms  involving  this  frequency  to  zero,  i.e., 

J2n  t  j2n  t 

ReaKKH^Cn^  )e  '  )G2(n,  ,  )e  }"0  (95) 

The  above  equation  describes  the  relationship  between  ,  H2,  K  and  G^ 
to  achieve  reduction  of  interference  at  the  frequency  2n^ .  Since  we 


have  chosen  as  the  frequency  for  compensation,  any  reduction  of  the 

interference  at  the  frequency  2Q^  would  be  incidental. 


Tf  reduction 


of  the  IM  distortion  at  both  2J2^  and  ^^2  required,  then  the 

problem  would  have  to  be  formulated  in  an  entirely  different  manner; 
perhaps,  the  center  frequency  of  the  compensator  along  with  its  bandwidth 
could  be  appropriately  chosen. 

Multiple  Frequency  Compensation 

We  have  stated  an  equation  which  describes  the  relationship  between 
the  sue  and  the  compensator  to  reduce  intermodulation  interference  at 
2i2^  ,  for  an  input  consisting  of  the  sum  of  two  unmodulated  tones.  Now 

let  us  compensate  the  IM  interference  for  the  case  of  an  input  which  has 
a  spectral  bandwidth  r.  Now  the  reduction  should  be  achieved  over  a 
band  over  which  IM  Interference  occurs: 

2SI^  ±  r,  instead  of  a  single  frequency  2ft^ . 

Quadratic  Mean-Square  IM  Criterion  and  a  Related  Cost  Function 

An  index  determining  the  amount  of  interference  reduction  is 
discussed  next.  It  is  natural  to  define  a  distortion  energy  criterion 
as 

00 

e  =  ~  /  |Y„(u)l^dw  (96) 

Ztt  2  ’ 

—  00 

The  actual  input  signal  x(t)  is  often  not  known,  except  that  it  is 
bandpass  and  real-valued  (as  is  also  the  linear  transfer  function 

H^(3)).  For  simplicity,  we  take  the  spectrum  of  the  input  to  be  uniform 
over  this  band,  so  that  Y2(a)^,(i)2)  “  Q2(w^.'»)2^’  since 

00 

=  — —  /  Q2(u)-4<^  ,  i|/^  )dij^^  (97) 

2it  -® 


46  - 


therefore  our  practical  distortion  criterion  becomes 


1  ^  p 

rr"  ^ 


where  A  is  the  frequency  range  of  interest. ( i.e.  The  region  in  which 

interference  is  to  be  minimized.)  It  will  be  seen  later  that  «  2a) 

0  c 

and  A  =  r  if  the  IM  interference  at  the  double  frequency  is  of  concern. 

The  cost  function  must  now  be  defined  over  the  appropriate  band  of 
frequencies.  Using  equation  (95)  the  intermodulation  at  the  (o)^  .a)^) 

frequency  pair  (in  the  associated-response  frequency-plane)  is 

Q2(a)^.a)2)  =  ^^(a)^  ,a)2)+H^  (w,  )H^  (a)2)G(a)^  .to^)  (99) 


r  r 

where  the  frequencies  and  a)^  each  vary  from  £1^  -  -  to  +  - 


Computation  of  the  Optimal  Parameters 

Parameter  optimization  applied  to  the  design  of  dynamic  systems 
Involves  the  selection  of  a  parameter  set  so  as  to  enable  the  system  to 
attain  the  beat  possible  performance.  We  have  selected  for  compensator 
design  the  nonlinear  least  squares  minimization  algorithm  by  Dennis  et 
al  [7].  This  algorithm,  in  the  large  residual  case,  is  more  reliable 
than  Gauss-Newton  or  Levenberg-Marquardt  methods  [8]  and  more  efficient 
than  the  secant  or  variable  metric  methods  algorithms  [9],  such  as  the 
Davidon-Fletcher-Powell  method  [10,11].  The  algorithm  amounts  to  a 
variation  on  Newton's  method  in  which  part  of  the  Hessian  matrix  is 
computed  exactly  and  part  is  approximated  by  a  quasi-Newton  updating 
method.  To  promote  convergence  from  poor  starting  guesses,  the 
algorithm  uses  a  model/trust-region  technique  along  with  an  adaptive 
choice  of  the  model  Hessian.  The  residual  vector  needed  by  NL2SN0  is 


2 

constructed  by  converting  the  interference  power  density  | ) | 

into  the  proper  form.  The  transformation  of  the  interference  density  in 
the  0)^  plane  into  the  real  world  frequency  axis  to  is  done  by 

searching  the  to^  .to^  plane  along  constant  sum  (or  difference)  frequencies 

of  the  argument,  and  summing  it.  This  is  expressed  by 

Qj^(to)  =  ;  Q^Cu^.to^)  dJ,  (100) 

to^  +(1)2=  u 

where  the  integral  is  a  line  integral  along  the  line  defined  by 


(0^ +0)2=0)  .  Note  that  dJl  and  do)  are  orthogonal  to  each  other.  The  10^10)2 

plane  along  with  the  incremental  line  element  dH  is  shown  in  Fig.  13. 

If  the  band  of  concern  were  the  baseband,  then  we  would  use  the  line 
defined  by  to^ -0)2=0)  for  the  above  integration. 

In  the  computer  program,  the  cost  function  is  simply  defined  as 


'PGM 


N  N 

i  \  u  \  l'=2^“li'“21<^l‘ 

i=-N  k=-N 


(101  ) 


where  o)^  ^  and  are  the  grid  points  in  the  (o)^  ,0)2)  frequency  plane 


r  r  r  r 

covering  the  square  2  ^  ^  ^^2~  2  ’  ^2*  2 


The  center  of  the  square  grid  is  (Q^,fi2)*  program  NL2SN0  uses 


r 


i,k 


A 


as  the  residual. 


|Q2(u)i..o,2k)|  (102) 

2 

The  dimensionality  of  the  residual  vector  is  (2N+1)  . 


Initial  Pareuneter  Selection 


In  most  optimization  methods  the  a  priori  values  affects  the 
convergence,  as  well  as  the  final  solution,  as  local  minima  may  exist 
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Fig.  13*  The  to^  ,a>2  plane. 


throughout  the  range  of  the  cost  function.  Therefore  a  method  of 
determining  suitable  initial  parameter  values  (even  in  a  semi-scientific 
way)  is  highly  desirable.  The  approach  we  use  is  to  compute  the 


parameter  values  needed  to  compensate  the  SUC  just  at  the  center 
frequency  (of  the  interference  band).  Then  equation  (99)  becomes 


0  . 


where  £2^  =£2^  »£22.  Rearranging,  we  obtain 


(103) 


’'h'J'h  ■  "hag  ® 


(104) 


Shown  in  Fig.  14  is  the  magnitude  of  Q2(£2^,£22)  in  the  >“2  Plane* 


■  ■  V  %'  V 


Consider  a  compensator  which  consists  of  a  linear  path  with  gain  K,  and 

a  quadratic  path  as  shown  in  Fig.  15.  The  transfer  functions  of  the 

-  49  - 


blocks  are 


G  (3)  = 
a 


s  +b 


s  *a. 


.  G^(3) 


3  +a„ 


From  this  the  NLTF  of  the  quadratic  path  becomes 


(s^  +a^  Xs^+s^+a^) 


(105) 


We  shall  call  this  a  Type  I  Quadratic  Compensator. 

Suppose  now  we  restrict  our  attention  to  a  single  unmodulated  carrier 

with  frequency  fl  and  the  interference  generated  at  the  frequency  2Q  . 

G  C 

Setting  3^  =3^  gives 


(106) 


To  compute  the  initial  parameter  vector,  we  equate  the  right-hand  side 
of  (1014)  and  (105).  We  will  now  compute  so  that  the  phase 

angle  of  the  compensator  is  equal  to  8„.  In  order  to  select  0„  so  that 

n  Ij 

it  is  the  additive  inverse  of  9„,  let 

n 

9q  =  9^  -  2ir,  if  9j^  >  0 


if  9^^  <  0 


Then,  equally  distribute  9^  among  the  compensator  poles  at  -a^  and 


"a2 

—  ,  and  zero  at  -b^  by  defining. 


Since  0„  is  a  negative  angle,  9_  is  positive. 

u  U 


The  initial  parameter  values  are  then  obtained  by  setting  (the  gain) 
bg  of  the  compensator  to  1  and  solving  expression  (106)  for  each 


parameter  value, 


(107) 


®D 

tan(^) 


®D 

tan(^) 


20 


tan(  Op) 


Using  these  computed  parameter  values  evaluate  the  compensator, 

J9u 


^MAG  ® 


setting, 


MAG 


MAG 


to  adjust  the  gain  parameter  from  its  previously  defined  value  of  1 


This  initial  parameter  method  provides  a  starting  location  for 
compensation  over  the  desired  bandwidth.  We  will  demonstrate  the  use  of 
this  result  with  an  example. 

Let 

0=5  rad/ sec 
c 

H  (0^)  =  0.01  -j 
I  c 

H^COc.^c)  *  -0.09654  -jO. 08098 
Then  from  (103)  the  compensator  response  should  be 
°2^”c’”c^  =  -0,0013  -j0,005 
which  gives 

-  -0.0013  ,  -  -0.005 

“mag”  =  -1.82516  rads 

and  since  0„  is  <  0  ,  0„  »  9- 
H  Hu 


/.v-V.Vl 


•  -•  «**'.• 

•  *.• 


0r,  »  0.60839 


-15.9269 


b  -  _ 2 _ 

‘^l  tan(  0.30M19) 

*1  ■  tan(  0.301119)  ”  '5-9269 

■  tan(  i?60839)  '  '“•'‘5* 

now  compute  the  using  these  parameters, 


2.2572 


“mag 

G - 

MAG 


The  initial  parameter  values  are  then 

£  -  (2.2572  -15.927  15.927  164. 165)^ 

EXAMPLES  ON  QUADRATIC  COMPENSATION  VIA  PROGRAM  VC0MP2 


We  now  give  two  examples  of  quadratic  interference  reduction.  Both 
examples  have  been  generated  by  use  of  the  computer  program  VC0MP2. 


Example  1 ; 

Consider  a  system  with  the  Linear  TF 
0.250025 


(s) 


s  +0.1S  +25.25 

which  has  a  bandpass  nature.  The  quadratic  nonlinearity  is  of  the  form 

100 


H2(3^ .s^) 


(s^  +2)(32  +2)(3^  +3^  +12) 


For  the  type  of  compensator  to  use  we  select  the  general  form  outlined 
in  the  previous  section. 


bo(3,+b,)(32  +  b^) 


(s^  +a^ 


The  input  is  a  narrowband  signal  with  a  center  frequency  of  5Mrad/3 


and  a  bandwidth  r,  of  2.0  Krad/s.  The  sum-squared  error  (i.e.,  sum 


squared  of  r(a)^  ,0)2)  of  the  system  with  no  compensation  over  a  5  x  5 

point  frequency  grid)  is  ep^=  3.9^727.  A  simpler  3x3  example  of  the 

frequency  grid  is  shown  in  Fig.  16.  Note  that  the  resulting  band  of 
interference  is  (2fi^“  r,  2n^+  r).  The  initial  parameter  vector  is 

p  =  [2.2572  -8.817  8.817  5.981  ]'^ 

After  optimization  the  compensator  parameter  values  are 
p  =  [^3.65  -129.9  2311.9  59.72] 

A  plot  of  the  sue  interference  function  20Log( I )  versus  the 

'  in  ' 


compensated  one  is  shown  in  Fig.  17.  The  system  quadratic  sum-squared 
error  is  EpQn|=  O.O817211,  a  vast  improvement  over  the  uncompensated  case. 


Fig.  16, 


3x3  point  frequency  grid  of  rCw^ .Wp) 


Example  2 : 

This  example  is  similar  to  Example  1,  except  the  amount  of  deviation  is 
increased  to  50.0  krad/s.  The  Initial  uncompensated  sum-square  error  of 
the  nonlinear  component  again  over  a  5  x  5  point  frequency  grid  for 
r((»).  ,0)-)  is  Ep^w-  8.301.  The  same  initial  parameter  values  as  in  the 


8.301.  The  same  initial  parameter  values  as  in  the 


last  example  are  used.  The  compensated  system  error  is  4.166.  The  final 
parameter  values  are 

£  -  [86.05  -4533  6480  190]”^ 

In  Fig.  18  a  plot  is  given  of  the  uncompensated  and  compensated 


•  v  v 
‘  '*'*'**] 


fci  •  •  ^ 

>  .'-•.'v  .V 
.■'■'.v'.V.v' 


Uncompensated 


p  •  5  ' 


10  dB 


Compensated 


4.997  4.999  5.001  5.003 


.  .N*. 


.  •  •  •  ^  .  *1 
-  I 

^  m  ^  ~  a  i 


Mr ad/s 


Fig.  17.  Compensated  versus  uncompensated  distortion. 


interference  function  20Log( | Q^„(a)) | )  for  this  example.  In  the  above 


examples  we  used  a  compensator  whose  blocks  were  single-pole  blocks.  A 
more  effective  reduction  of  the  IM  interference  could  be  achieved  with  a 
more  complex  conpensator  whose  block  TFs  are  second-order.  Its  NLTF  is 
given  by. 


G2(s^  .s^) 


_ t>o(3,>b^)(3g^b3) _ 

(s^+a^s^+a^ )(S2+a232+a^ ) ( ( s^ +S2)^+a2( 3^ +S2)+a^ ) 


(108) 
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Fig.  18.  Compensation  results  of  the  second  example 

This  form  of  compensator  will  be  called  Type  II  Quadratic  Compensator. 

In  this  section  we  have  explored  the  compensation  of  interference 
produced  by  a  quadratic  nonlinearity.  As  stated  earlier,  the 
Interference  generated  by  quadratic  nonlinearities  does  not  fall  in  the 
system's  passband.  However,  the  Intermodulation  terms  generated  by 
cubic  nonlinearities,  can  fall  within  the  main  passband  of  the  input  to 
the  sue.  Compensation  of  third-order  nonlinearities  is  considered  next. 

COMPEHSATION  OF  CUBIC  NONLIMEARITIES 

In  this  section  we  present  the  design  of  the  third-order  nonlinear 
compensator  using  the  configuration  as  shown  in  Fig.  19.  The  SUC  here 
is  shown  in  greater  detail  in  Fig.  20.  Note  that  its  cubic  nonlinearity 
may  be  represented  mathematically  by  the  NLTF 


H3^®1,®2'33)  “  "a^3,)H^(32)H^(33)H^(3^. 3^-33) 


(109) 


i 


Consider  as  an  input  to  the  SUC,  a  signal  consisting  of  the  sum  of  three 
sinusoids, 

x(t)  =  cosset  +V^cos£j^t 

the  output  of  the  SUC  is  the  sum  of  the  outputs  pf  the  two  paths;  the 
response  of  the  linear  block  is 

w^(t)  =  vjH^(n^)|cos((a^)t 


I 


+V2|H^(n2)|co3((n2)t 


+V^|H^ (n^) lcos((ft^)t 


(110) 


SYSTEM  UfJOER  COMPENSATION  COMPENSATOR 


i 


Fig.  19.  Post-compensation  of  a  communications  system 
with  cubic  nonlinearity. 


The  response  of  the  cubic  block  using  expression  (29)  is 
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HH,(3«2)ho3((3n2)t 

HH^(2ft^>fl3)lcos({2a^^a3)t  .26^^-  e,o^"<t>2n, 

. (H^(2fl,  -«3)  1 003((2a^  -Q^n  ~\^2a^ 

-r,23lQ,23l[|Hd(«rV3^l^°®^^“rV“3^^  "\'%'VS"«2'“ 


>lH^(fi^-fl2-fl3)lco3((a^-ft2>a3)t  "e£j^-6n^"en^"l>n^-ft2-a 


*ri33|Qi33lK(nr2«3^l®°®^^“r^”3^^ 

+2|H^(a^)lcos((a^)t  +2ejj  ) 

.  I  H^(  a,  *2113)  I  oo3(  ( a,  .aOj)  t  .  ‘a,  .2a’  ^ 

*''223l«223l(=l"a<”3>l“*<<”3’‘ 

.|H  (2a  •a.)|aos((2a2.a3)t  *2eo*»a  *♦28  .a 


2  “3  ‘-"2  3 


-r233lQ233l[|H^(V2«3^l"°®^^V2«3^" 

+2|H^(a2)|co3((a2)t  +26^^  ) 

^IH  (a  -2£j  ) I C03((a2 -2a  )t  -2e  2a,^^ 


2  3  2  ”3 


^^333'®333  ^^^33^  3  a3 
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+  lH^(3fi2)|co3((3fi3)t  +30JJ  )]}  (111) 

The  spectra  of  and  w^,  the  output  of  the  SUC  for  an  input  consisting 

of  the  sum  of  three  sinusoids  is  shown  in  Fig.  21.  The  in-band 
distortion  falls  at  the  following  frequencies:  ^2 

.  In  the  ensuing  discussion  we  have  chosen  the  distortion  most 

adverse  to  our  communication  system  to  be  at  a  band  of  frequencies 
around 


Post-Compensator  Representation 

The  cubic  compensator  structure  shown  in  Fig.  22  appeared  promising 
and  was  tried  first  in  our  studies.  The  parallel  paths  of  this 
compensator  are  a  gain  constant  for  the  linear  TF  and  a  third-order 
nonlinear  path.  Only  the  nonlinear  response  of  the  SUC  will  be 
compensated;  the  linear  response  will  be  left  essentially  unaltered. 


Fig.  21 .  Spectral  plot  of  a  third-order  SUC  output 


Fig,  22.  Cubic  nonlinearity  compensator. 

The  output  frequencies  generated  by  the  channel  NLTF  are  out-of-band  and 
in-band  interference  terms.  Our  choice  of  an  interference  frequency  to 
be  reduced  is  in-band  term.  Having  chosen  the 

interference  frequency  of  interest,  the  compensator  design  will  now  be 
discussed.  Note  that  at  the  moment  we  are  dealing  with  the  case  of  an 
input  consisting  of  three  sinusoids. 

The  output  of  the  compensator  is 

y  =  yT>y3 

where  the  output  from  the  linear  section  of  the  compensator  is 

y^  =  KG^ [w^ ]+KG^ [w^]  +  •••  (112) 

^1  =  ^I.L  ^  ^I.N 

The  compensators  cubic  output  is 

y^  =  G^[w^  ]+G^[w^]+  3G^{w^  ,w^  ,w^ }+  3G^(w^  ,w^,w^)  +  •••  (113) 

Since  our  analysis  is  restricted,  for  the  moment,  to  the  third-order 
response,  we  will  omit  all  terms  shown  underscored.  Now,  substituting 
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for  and  from  (110)  and  (111)  and  retaining  using  only  the 
distortion  terms  at  the  frequency  of  interest  we  have 


=  +(j)^  )]  (iu) 


‘IM 


=  G3[w^] 


‘IM 


=  ev^v^v^lH^ (n^ )H^ (-fl2)H^(a3)G3(fti  .-n2.«3)| 


•003(  (i!,  -a^.a^) 


(115) 


where  we  have  let  be  a  scalar  gain  K. 

From  earlier  discussion,  we  know  that  removal  of  the  third-order 
output  effects  requires 

°1»1  =  ' 

G1H3.G3H1  .  0 

or,  in  terms  of  the  composite  system  TFs, 

Ql  =  I 

Q3  =  0 

Since  we  stated  that  our  goal  is  to  reduce  the  IM  interference  at 
we  do  so  by  setting  the  sum  of  the  terms  involving  this 

frequency  to  zero,  i.e.. 


6KV,V3V3lH3(a,,-0j,i.3)|co3((!l,-V‘)3)t  *  ♦a,-aj.a3> 

♦6v,  VjV Jh,  (a,  )H,  (-QjjH,  (83)0  (a,  .-Oj.ii  )  | 


•O03((a,-ii2*03)t  •ea/»-a/“«3*\ -0^.83’  '  “ 


(116) 


The  above  equation  describes  the  relationship  between  H. ,  H.,  K,  and  G 
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to  achieve  reduction  of  interference  at  the  frequency  Q^u  -fl, -a„+a-. 

iPl  1^3 


Multiple  Frequency  Conpensatlon 

We  have  stated  an  equation  which  describes  the  relationship  between 
the  sue  and  the  compensator  to  reduce  intermodulation  interference  at 
a^-a^+a^,  for  an  input  consisting  of  the  sum  of  three  unmodulated  tones. 

Now  let  us  consider  the  IM  interference  for  the  case  of  an  input  which 
has  a  spectral  distribution  over  a  bandwidth  r.  An  example  spectral 
distribution  of  the  channel  NLTF  response  is  given  in  Figure  23. 

Now  the  reduction  should  be  achieved  over  a  band  over  which  the  IM 

3r 

interference  occurs:  (a^-a^+a^)*  ^  instead  of  merely  the  frequency  a^- 

a2+a2. 

Cubic  Distortion  Criterion  and  a  Related  Cost  Function 

A  cost  function  describing  the  amount  of  interference  is  formulated 
next.  It  is  natural  to  define  a  distortion  energy  criterion  as 


1_ 

2it 


/  |Y^(w)|  dto 


(117) 


The  actual  input  signal  x(t)  is  often  not  known,  except  that  it  is 
bandpass  (as  is  also  the  linear  transfer  function  H^(s)).  For 

simplicity,  we  take  the  spectrum  of  the  input  to  be  uniform  over  this 
band,  so  that  )  =•  ) .  Now  if  we  define 


Q2(  id) 


(2tt)' 


/  /  Q^(u)-i|(,  ,  d'1'2 


(118; 
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OUTPUT 


Fig.  23.  Spectral  plot  of  a  third-order  SUC  output  from  an  input 
consisting  of  a  carrier  with  bandwidth  r. 


then  our  practical  distortion  criterion  becomes 


1  ^  p 


(119) 


where  i*)q±A  is  the  frequency  range  of  interest. (i.e.  The  region  in  which 

interference  is  to  be  minimized.)  It  will  be  seen  later  that  u)„  »  n 

0  c 

and  A  =3r  where  the  IM  interference  at  SJ-u  “R, is  addressed. 

IM  123 

The  cost  function  must  now  be  defined  over  the  appropriate  band  of 
frequencies.  Using  equation  (118)  the  intermodulation  at  the  (u^ 


frequency  coordinates  (in  the  associated-response  frequency-plane)  is 


( u)^  , I (i»2 )  =  KH2(a)^ 


(120) 


r  r 

where  the  frequencies  to^  ,  oj^  and  each  vary  from  (n^-  -  )  to  (n^+  -  ). 

The  residual  vector  needed  by  NL2SN0  [14]  is  constructed  by  converting 

the  interference  density  1 02(01^  I  into  the  proper  form.  The 

transformation  of  the  in-band  interference  power  density  in  ( , lo^ 1 ) 

into  the  real  world  frequency  m  is  accomplished  by  integrating  the 
values  of  that  lie  on  an  interference  surface  Ojh^Cid)  defined  by 

(1)^  -  (4)2  +  (0^  =  OJ 

This  is  expressed  by 

Qim(w)  =  If  , W2 . Wj )  doj^  (121) 

u)  ”0)^  -W2'''<D2 

where  the  integral  is  a  surface  integral  along  the  plane  defined  by 

Otu(u)  i’on  u), -a)„+ 0)^-0)  .  Note  that  do^w  and  du  are  orthogonal  to  each 
In  123  in 

Other. 

In  the  computer  program,  the  cost  function  is  simply  defined  as 


N  N  N 


'PGM  -III 

i--N  j=-N  k=-N  ^  ^  ^ 


(122) 


where  w  ,  w  and  u-.  are  the  grid  points  in  the  (w  ,u)_,w,)  frequency 

1  1  <1 J  j  K  1  t  j 

r  r  r  r  r 

space  covering  the  cube  [n^-  -  ,  2  ^  x  -  .  "^2+  ^  ]  x  -  , 

r 

-  ].  The  center  of  the  cubic  grid  is  program  NL2SNC 

uses 


^.j.k  =  IS^'^1i'“2j'“3k^l 


i 


(123) 


as  the  residual.  The  dimensionality  of  the  residual  vector  is  (2N+1 )  . 

Shown  in  Fig.  24  is  a  plot  of  the  r.  .  grid.  Fig.  25  shows  the 

i  »  J  » 

location  of  an  example  cube  centered  at  ( ) . 

Fig.  26  shows  one  of  the  surfaces  the  three-dimensional 

IM  G 

interference  space  for  this  plane  maps  into  the  real  world 

frequency 


Delay  Block  In  the  Compensator  Structure 


Let  us  define  the  group  delay  of  the  contributions  of  the  third- 
order  intermodulation  as 

T  (u)  =  “  Phase(Q,„(u)) }  (124) 

oto 
IM 


We  first  considered  the  cubic  compensator  previously  shown  in  Fig. 
21  with  nonlinear  path  NLTF  given  as 


G^(s^  ,3^,3^)  = 


b^Cs^+b^  ^ 


(125) 


(s^  +a^  ■^S2^®3'^°2^ 


which  we  will  call  the  Type  I  compensator .( A  fairly  complete  family 

listing  of  the  useful  cubic  compensator  structures  is  given  in  Table  2.) 

When  the  group  delay  of  the  third-order  nonlinearity  of  the  SUC  was 

small,  the  optimized  Type  I  compensator  turns  out  to  be  stable  (i.e., 

the  block  poles  of  the  compensator  turned  out  to  be  in  the  left-half  s- 

plane).  However,  large  values  of  SUC  group  delay  t  render  the  Type  I 

IM 


compensator  inadequate;  for  such  cases  the  program  VC0MP3  yields  a)  only 

a  small  reduction  in  or,  b)  unstable  block  poles  for  the 

PCjM 


compensator,  or  both. 


Type  I 


K 


(s  +  b^) 
(3  +  ) 


(s  +  a^) 


Type  Ip  The  same  aa  the  Type  I;  except  that  =  Ke 


-ST 


Type  Ila 


bQS  Tb^s.  b3 


(s  +  ion  )(3  +  o.in  ) 
c  c 


(3  +  a^) 


Type  11 Sp 


The  same  as  the  Type  Ila;  except  that  »  Ke 


Type  lib  K 


bpS  ^  b,3T  b3 


s  T  a^sT  ap 


1 

(s  +  a3) 


Type  Ilbp  The  same  as  the  Type  lib;  except  that  G^  =  Ke 

Note:  Structures  I,  Ip,  Ilb,  and  Ilbp  are  not 

guaranteed  to  be  stable.  However,  as  discussed 
in  subsection  3.4,  they  can  be  reparametrized  so 
that  in  the  space  of  the  new  parameters  these 
forms  are  guaranteed  to  be  stable  for  every 


Fig.  an,  r  grid 

1 1  J  # 

A  modified  compensator  structure  which  can  effectively  reduce  the 
IM  interference  is  proposed  next.  This  modified  structure  contains  a 
delay  block  in  the  linear  path  of  the  compensator,  as  shown  in  Fig.  27. 
The  compensator  NLTF  remains  the  same,  thus  the  designation  of  this 

compensator  is  Type  The  Incorporation  of  the  delay  shifts  the  phase 

of  w^  (sue  IM  interference)  so  that  the  group  delay  of  y^  (the 

compensator  cubic  response  to  the  input  of  w^ )  duplicates  the  group 

delay  of  y^  delayed  SUC  IM  interference).  This  altering  of  the 

group  delay  facilitates  constructing  a  stable  so  that  y^  is  the 

additive  inverse  of  y^  The  delayed  version  of  w^  (i.e.,  y^  j^),  has  a 

phase  characteristic  that  can  be  adequately  cancelled  by  the  response  y^ 
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Fig.  26.  r.  grid  with  one  interference  surface  shown, 
i  I J  » ^ 


unstable  block  poles.  The  output  of  the  compensator  is 


y  =  yi,L  "yjM  * 


(126) 


where  is  the  linear  response,  and  is  the  IM  interference  term, 


and  the  underscored  terms  are  of  order  higher  than  three,  hence  omitted. 


Note  that 


yi,L= 


where  is  the  delay  operator,  D^[w^]  =  w^(t-T) 


(12T) 


Clearly,  the  IM  interference  signal  is 


y^(t)  =  KD^[H^[x(t)]]*  G^Cw^(t)]. 


(128) 


The  effect  of  the  group  delay  manifests  into  the  composite  system  NLTF 


in  the  following  manner: 


-(s^  +S2+S  )t 

(s^ ,S2,s^)  =  Ke  ^  H^(3^ ,S2,s^) 


+  H(s^  )H(32)H(3^)G^(s^ ,32,3^) 


(129) 
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Fig.  27.  Addition  of  a  phase  delay  to  the  cubic  compensator 


The  Type  [._  cubic  compensator  provided  an  improvement  over  the  Type 
D 

I;  the  computer  program  VC0MP3  returns  unstable  block  TFs  only 
infrequently.  Still  when  compensating  an  SUC  with  a  large  group  delay 
it  was  occasionally  unstable.  This  led  to  the  development  of  yet 
another  type  of  cubic  compensator  in  which  the  poles  of  the  block  TFs 
are  restricted  to  only  stable  values.  In  this  compensator  the  form  of 
the  block  TFs  is  the  same  as  in  the  Type  II;  however  the  denominator 

parameters  are  now  frozen  as 
'> 

a„  =  «■ 

0  c 

a  =  o.in  +  ion 

1  c  c 

30  that  the  the  compensator  poles  are  at  lOn^and  O.in^.  Thus, 


G^(s)  . 


.^q _ 

(s+  (o.i)n^)  (s*  (iO)n,J 


,  G^(s)  =  1 
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(122) 
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This  structure  is  called  Type  Ila  cubic  compensator.  We  now  discuss  a 
method  to  allow  the  compensator  block  poles  to  vary  while  simultaneously 
ensuring  that  they  stay  in  the  left  half-plane. 

Re-parmetrizatlon  of  the  Compensator  Block  TFs 

To  completely  circumvent  the  possibility  of  obtaining  unstable 
poles  in  the  compensator  block  TFs,  we  re-parametrize  the  G  and  G 

3  G 

block  TFs  in  G^.  For  brevity  we  will  only  consider  Type 
structure.  The  TF  structure  for  G  (s)  and  G  (s)  is  expressed  as: 

3  Q 


:  v’y-^ 

» v"’- 

If. 


The  new  parameter  vector  is 


wherein  only  P^.P^  affect  the  poles  of  G^,  and  Pg  affects  the  pole  of 
G^.  Every  point  of  this  parameter  space  maps  into  a  stable  point  in  the 
(pQ ,b^  ,a^  .Sq , c, t)  space.  Indeed  it  can  be  seen  that  the  damping 

ratio  of  G  (s)  is  guaranteed  to  be  greater  than  or  equal  to  i;.  This 


s^  and  lie  to  the  left  of  the  contour  T  no  matter  what  the  values  of 


-ji^ 


P  “  ^  P-j  » p2  *  P^  » P^l  >  P^  »  Pg  » Py  ^ 


form  will  be  called  the  Type  The  poles  of  the  Type  Hb^ 

compensator  are  shown  in  Fig.  28.  Notice  that  the  compensator  poles  at 


the  parameters  p.,  and  p^  may  be. 

4  b 


Fig.  28.  Pole  plot  for  modified  Type  Ilb  cubic  compensator  designs. 


EXAMPLES  OF  INTERMODULATION  REDUCTION  USING 
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THE  PROGRAM  VC0MP3 


Results  of  example  runs  on  the  program  VC0MP3  are  presented,  which 
confirm  the  effectiveness  of  our  technique  in  reduction  of  IM 
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interference  in  a  specified  band. 


Since  most  real  world  communication  channels  exhibit  band-pass 
behavior,  the  model  used  for  the  channel  is  a  Chebyshev  bandpass  filter 
of  second,  fourth  and  sixth-degree  (passband  ripple  is  0.5  dB).  The 
nonlinearity  is  third-order  and  produces  a  signal  approximately  -20  dB 
below  the  linear  response.  In  the  computer  examples  we  vary  the  cutoff 
rate  of  the  channel  (degree),  the  NLTF  of  the  channel.  The  compensator 
used  is  of  Type  11^^  form  shown  earlier. 


Compensation  Excunples  by  the  Computer  Program  VC0MP3 

We  now  give  "everal  numerical  examples  of  cubic  interference  reduction 
by  the  computer  program  VC0MP3. 


Example  1 

Consider  a  system  with  the  transfer  function. 


(s)  = 


0.571997(10^)3 


s^+  0.571997(10^)s  +24.9916(10^^) 


which  is  a  second-degree  Chebyshev  with  a  passband  4.9  Mrad/s  to  5.1 
Mrad/s.  The  system  has  a  cubic  nonlinearity  of  the  form. 

H^(s^  ,82,3^) 


4.6416 

(s^+  2)(s2+  2)(s^+  2)(3^+S2+s^+  12) 


This  NLTF  has  a  three  dimensional  impulse  response  resembling  a  lowpass 
filter.  The  type  of  compensator  used  is  the  stable  cubic  with  linear 
delay,  outlined  in  the  previous  section. 


03(3^,32,33,2) 


(PjS^ 


P,s,* 


(s  +  0.1a)„)(3  +  10a)_)(3.+  O.lo)  )(s„+  10a)_) 


^P2VPlV 


Pq^ 


(3^+  0.1(u^)(3^+  lOtu^) 


with  the  linear  eection  of  the  compeneator  defined  as 

-(3  ^-3  +3  )p 
G(s)  =  e  ^  ^  i  ^ 

The  input  is  a  narrowband  signal  with  a  center  frequency  of  5  Mrad/s 

and  a  bandwidth  r  of  200  Krad/s.  The  sum-squared  error  of  the  system 
with  no  compensation  is  tp^  =  0.30367*48.  The  initial  parameter  vector 

is 

£  -  [-*(. 6*116  1  .587  10.33  0.0]^ 

After  optimization  the  compensator  parameter  values  are 
P  =  [-839. 3  -32.18  -33.21  1.265] 

The  poles  of  the  compensator  were  fixed,  only  the  gain  and  the  location 
of  the  zeros  were  varied.  We  therefore  know  that  this  compensator  is 
stable.  The  compensated  system  Cp^  is  0.0003306  a  vast  improvement 

over  the  uncompensated  case.  A  plot  of  the  compensated  nonlinear 
interference  verses  the  uncompensated  is  given  in  Fig.  29- 
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A. 8  4.9  5.0  5.1  5.2 

Mrad/s 


Fig.  29.  Compensated  versus  the  uncompensated  Interference. 

Example  2 

This  example  is  similar  to  Example  1.  The  linear  channel  TF  is 
unchanged;  however,  the  blocks  of  the  NLTF  are  all  second-degree 
bandpass  Chebyshev.  Thus  =H^  both  have  a  somewhat  narrow  bandpass 

characteristic.  The  sum-squared  error  of  the  system  with  no 
compensation  GpQ^  is  105.3^013.  The  initial  parameter  vector  is 

p  =  [-4.6416  0.0  0.0  0.0]^ 

After  optimization  the  compensator  parameter  values  are 
p  =  [1647.3  -85.75  66.65  1.579] 

The  compensated  system  is  0.0467988,  again  a  vast  improvement  over 

the  uncompensated  case.  A  comparison  of  the  compensated  nonlinear 
interference  and  the  uncompensated  one  is  given  in  Fig.  30 


Mrad/s 


Fig.  30.  Compensated  versus  the  uncompensated  interference. 


Example  3 

This  example  has  a  linear  TF  with  a  further  reduced  bandwidth.  Th( 
linear  characteristic  of  the  channel  is  modeled  by  a  fourth-degree 
Chebyshev  with  a  passband  of  4.9  Mrad/s  to  5.1  Mrad/s. 

H^(3)  - 

_ 0.057144073^ _ 

[s**+  0.284847323^+  50.04374633^+  7.118792693  +  624.580478] 

The  sum-squared  error  of  the  system  without  compensation  is  epQ„] 
=59.27650.  The  initial  parameter  vector  is 
p  =  [-25.0  0.3  -1 .0  0.0]"^ 

The  optimally  compensated  system  is  0.06369799.  The  final 

parameter  values  are 

p  =  [-3925  -22.78  -157.0  0.002317]^ 


A  plot  of  the  compensated  nonlinear  interference  as  well  as  the 
uncompensated  one  is  given  in  Fig.  31 .  Notice  that  20dB  of  interference 
reduction  has  been  achieved  with  only  a  second-order  compensator. 


4.0  5.0  6.0 

Mrad/s 


Fig.  31 .  Compensated  versus  the  uncompensated  interference. 

Example  4 

This  example  has  a  linear  TF  with  a  very  narrow  bandwidth.  The 
linear  characteristic  of  the  channel  is  modeled  by  a  sixth-degree 
Chebyshev  with  a  passband  of  4.9Mrad/s  to  5.1 Mrad/s. 

H^(3)  = 

_ 0.0057083^ _ 

3^+  0.2503383^+  75.036103*^+  12. 518433^+  1875.2723^+  156.3556s  +  15609 

The  initial  parameter  vector  is 

£  =  [-38.07  0.2547  1.522  0]^ 

The  uncompensated  system  is  68.5240325.  The  final  parameter  values 

are 

£  =  [11390  37.12  455.8  2.2]^ 

The  system  MSE  is  21 .5367901  representing  a  narrow  improvement  over  the 


CHAPTER  V 


CONCLUSIONS 

The  performance  of  a  communication  system  can  suffer  degradation 
due  to  the  intermodulation  (IM)  interference  created  by  the  inherent 
nonlinearities.  It  has  been  shown  here  that  such  IM  interference  can 
be  effectively  reduced  for  the  frequency  band(s)  of  interference  by  a 
stable  post-compensator  appended  to  the  System  Under  Compensation  (SUC). 
A  complete  methodology  for  designing  such  a  compensator  has  been 
developed.  It  is  based  upon  a  mean-square  IM  criterion  and  was 
successfully  implemented  in  a  FORTRAN  program  VC0MP3.  This  program  uses 
a  highly  powerful  software  package  for  optimization  of  mean-square 
nonlinear  functions,  thereby  ensuring  convergence  to  the  global  minimum 
in  almost  all  practical  cases. 

The  key  strengths  of  the  new  approach  are  (a)  broadband 
compensation  (over  the  band  of  interest)  in  contrast  with  the  single 
frequency  compensation  attempted  in  the  past,  and  (b)  guaranteed 
stability  of  the  compensator.  The  latter  was  achieved  in  a  highly 
innovative  way.  The  stable  region  of  the  space  of  the  parameters  (of 
block  transfer  function  of  the  compensator)  was  mapped  onto  a  new 
parameter  space.  In  fact  the  stable  region  was  the  interior  of  the  s- 
plane  to  the  left  of  the  constant  damping  ratio  lines,  with  a  damping 
ratio  5  specified  by  the  designer.  Thus  for  the  reparameterized  problem 
the  program  VCOMP3  is  unconstrained  to  choose  any  point  in  the  new 


f. 


K 


Cr 

u 


space.  The  final  design  is  not  only  stable,  the  block  transfer 
functions  have  a  damping  ratio  greater  than  5. 

The  compensator  structures  proposed,  and  used,  have  blocks  with 
pole-zero  transfer  functions,  hence  only  low  order  blocks  are  required. 
Indeed  as  shown  in  the  studies  presented,  only  first  or  second-order 
blocks  were  required.  The  concept  of  group  delay  for  the  interference 
component  of  the  SUC  was  defined,  which  was  followed  by  an  improved 
compensator  design.  The  key  advantage  realized  is  that  the 
dimensionality  of  the  parameter  vector  was  small,  typically  3  to  9, 
depending  on  the  particular  structure.  Several  case  examples  were 
presented  which  demonstrate  that  reductions  of  15  to  50  dB  in  the  IM 
interference  can  be  achieved. 

For  completeness  of  this  volume,  a  fairly  comprehensive  discussion 
of  Volterra  systems  was  included.  This  research  clearly  shows  that 
advanced  theory  can  be  utilized  to  yield  a  practical  tool  for 
improving  the  performance  of  expensive  communication  systems. 

Future  work  can  entail  a)  simultaneous  IM  interference  in  several 
bands,  and  b)  simultaneous  IM  interference  reduction  for  multiple 
receivers. 
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APPENDIX 


SIMPLE  LINEAR  COMPENSATION  EXAMPLE 


In  this  section  we  offer  a  tutorial  example  in  altering  linear 
system  response.  This  material  is  essential  to  the  compensating  of 

and  it  provides  us  with  a  fundamental  compensator  structure  to  reduce 
the  effects  of  the  higher  order  terms. 

These  examples  deal  with  the  most  rudimentary  adjustment  of  linear 
distortion  in  a  communication  channel.  In  both  examples  we  take  the 
general  form  of  the  compensator  to  be 


of  interest  is 
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The  system  frequency  characteristics  are  given  in  Table  3.  As  a  measure 
of  the  effect  of  the  compensator  on  the  phase  characteristic,  a  linear 
approximation  to  the  phase  of  each  example  will  be  given. 

Linear  Approximation  to  the  slope 
m»a)  +c 

m  =  -0.889-10~^ 
c  =  -i4.289 

Maximum  phase  deviation  from  a  linear  approximation 
3.367° 

TABLE  3 

THE  UNCOMPEHSATED  SYSTEM  FREQUENCY  CHARACTERISTICS 


1^6 

(0  X  10 

Magnitude 

Phase (degrees) 

0.90 

0.177 

89.87 

0.95 

0.244 

39.30 

1.00 

0.249 

-4.28 

1 .05 

0.239 

-46.81 

1.10 

0.176 

-89.91 

It  is  desired  to  provide 

a  second-order 

post  compensator  to  flatten 

magnitude  characteristic 

over  this  band. 

The  compensator  structure 

chosen  because  of  its  wide  bandwidth  coverage  is 


G^(s) 
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K  is  chosen  to  provide  unit  gain  at  the  center  frequency  Thus  is 

the  only  free  design  parameter.  By  experimentation  a  satisfactory  value 
was  found  to  be  b^  =  200.0,  so  that  K  =  50.6^.  The  corresponding 

magnitude  and  phase  of  the  compensator  and  the  corrected  system  are 
given  in  Table  i*.  and  Table  3. 

TABLE  4 

THE  COMPENSATED  SYSTEM  FREQUENCY  CHARACTERISTICS 

The  Compensator  Response 


0)  X  10^ 

Magnitude 

Phase(degrees) 

0,90 

5.686 

-43.92 

0,95 

4.464 

-25.46 

1  ,00 

4.01 1 

0.00 

1  ,05 

4.423 

24.37 

1.10 

5.420 

41.19 

The  final  system  response  is  shown  in  Table  5,  overall  response  is  a 
flat  pass- band. 

Linear  Approximation  to  the  slope 
m*u)  ■'■c 

m  =  -0.433.10"^ 
c  =  -4.289 

Maximum  phase  deviation  from  a  linear  approximation 
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TABLE  5 

COMPEMSATED  SYSTEM  FREQUENCY  CHARACTERISTICS 

O 

X 

The  Final  System  Response 

Magnitude 

Phase(degrees) 

0.90 

1 .008 

45.92 

kS 

0.95 

1 .092 

13.83 

r\'\' 

1 .00 

1 .000 

-4,28 

1  .05 

1  .058 

-22.43 

1.10 

0.994 

-48.72 

1 

The  compensator  parameter  values  which  provided  the  desired  results  are 
listed  below: 

ao  =  1.0 

a^  -  10.1 
b^  =  1930.0 
b^  =  J405.0 
b^  =  193 
K  =  50.6^1 

Several  observations  can  be  made  concerning  this  example;  first,  the 
compensated  magnitude  characteristic  is  quite  flat.  Indeed,  it  is  found 
that  the  maximum  deviation  from  OdB  is  only  0.2dB.  Second,  note  that 
the  maximum  deviation  of  the  uncompensated  phase  from  a  linear 
approximation  given  by 

i>„(u))  =  -0.89066(10  ^)w  -4.28  (in  degrees) 

H 


is  3**^3  •  For  the  compensated  system,  the  maximum  phase  deviation  from 
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the  linear  approximation 

-  (0.015326)01 

is  0.878°. 


Example  2 


We  consider  the  same  communication  system  as  before  however,  here 


the  objective  is  to  seek  a  compensator  to  reduce  the  3dB  band-width  from 
200  Krad/s  to  100  Krad/s.  A  compensator  structure  is  chosen  that  will 


narrow  the  pass- band  it  is 


G^(3) 


(s  +(0.1)o)q) 

2  2 
s  +a^  s  +oiq 


As  before  K  is  chosen  to  provide  unit  gain  at  the  center  frequency 

The  magnitude  and  phase  of  the  compensator  is  given  in  Table  6.  The 
magnitude  and  phase  of  the  compensated  system  is  given  in  Table  7. 


TABLE  6 

THE  COMPENSATED  SYSTEM  FREQUENCY  CHARACTERISTICS 
The  Compensator  Response 


U)  X  10^ 

Magnitude 

Phase (degrees) 

0.90 

1 .813 

61  .92 

0.95 

2.894 

43.22 

1 .00 

4.010 

0.00 

1 .05 

2.950 

-42.09 

1.10 

1.940 

-59.87 

The  system  overall  is  a  3dB  narrower  pass-band. 
Linear  Approximation  to  the  slope 


W  ’irt  •  ’ 


in*(i)  *c 

m  -  -1.600-10“^ 
c  -  -4.289 

Maximum  phase  deviation  from  a  linear  approximation 


14.49 


TABLE  7 

THE  COMPENSATED  SYSTEM  FREQUENCY  CHARACTERISTICS 
The  Final  System  Response 


1  X  10^ 

Magnitude 

Phase( degrees) 

0.90 

0.321 

151.80 

0,95 

0.708 

82.53 

1 .00 

1.000 

-4.28 

1 .05 

0.705 

-88.90 

1.10 

0.343 

-152.01 

The  compensator  parameter  values  which  provided  the  desired  results  are 
listed  below: 

ao  »  1.0 

a^  -  0.107 
b^  -  0,0421 
b^  -  0.425 
b^  -  0.421 
K  -  50.31 

Several  observations  can  be  made  concerning  this  second  example;  first, 
the  compensated  magnitude  characteristic  reaches  the  -3dB  points  in  half 
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the  bandwidth  required  in  the  original  system.  Second,  note  that  the 
maximum  deviation  of  the  uncompensated  phase  from  a  linear  approximation 
given  by 


4>p^(a)) 


-1 .600(10”^)u  -14.28 


(in  degrees) 


is  1  <4. 14°. 
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